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Abstract Measurements have been made of the entrainment of water droplets in an air-stream 


when air and water flow concurrently up a vertical glass tube of 0-5in. diameter. It was found 


that the mass flow of entrained liquid varied with the air rate to the power 2-6 and the liquid 


rate to the power 1-5 when annular flow conditions prevailed. The effect of entrainment on the 


calculation of the thickness of the annular liquid film by the method given in Part I is discussed. 


Réesumeé— Les autuers ont fait des mesures sur lentrainement de gouttelettes d'eau par un 


écoulement dair quand leau et lair se déplacent dans le méme sens et vers le haut dans un tube 


de verre vertical de 0.5 in. de diamétre. Ils ont trouve que Vécoulement massique du liquide 


entrainé varie avec la vitesse de lair a la puissance 2,6 ct la vitesse du liquide & la puissance 1.5 


quand les conditions d’écoulement annulaire prevalent. L effet de lentrainment sur le calcul 


de Tépaisseur du film liquide annulaire par la méthode donne partie I est discute 


Zusammenfassung-—Es wurden Messungen durchgefiihrt tiber die Mitnahme von Wassertropfen 


durch einen Luftstrom, wobei Luft und Wasser im Gleichstrom aufwiirts durch ein senkrechtes 


Glasrohr mit 0.5 Durchmesser strOmten. Dabei wurde gefunden, dass der Mengendurchsatz der 


eingegebenen Flissigkeit sich mit der 2.6ten Potenz der Luftgeschwindigkeit und dass sich die 
Fliissivkeitsgeschwindigkeit mit der 1,5ten Potenz indert, wenn die Bedingungen einer Ring- 
strOmung herrschen, Der Einfluss der Mitnahme auf die Berechnung der Dicke des Ring-Fliissig- 


keits-Films mit Hilfe der in Teil I angegebenen Methode wird diskutiert. 


INTRODUCTION from the liquid flowing in the annular film, the 


il] on the pressure eradient arrangement 1, Fig. 1 was used. A short brass 


IN A previous paper 


both tube T fits closely round the glass tube GT in 


and hold-up in a vertical tube carrying 


which air and water are flowing cocurrently 


liquid and gas, attention was drawn to the fact 
upwards, Tube T carries a copper dome D fitted 


that a proportion of the liquid travels as droplets 
with a circular knife edge KE 0-8 times the inside 


within the gas phase. Although this phenomenon 


diameter of the glass tube and lin. lone. The 


has been observed by many workers, little quan- 


central core of fluid leaves uninterrupted through 


titative work has been reported, This paper 


describes a method of measuring quantitatively KE and the rest including the liquid film leaves 


through E and returns to the feed tank. The 


the local mass flow density of the liquid drops 
function of the flow probe PL is used to sample the gas core ; it is 


and their velocity as a 


connected via tube MC to a differential manometer 


variables, in a glass tube of approximately 


0-5in. diameter. The influence of this liquid en- Which indicates the impact pressure on the probe 


trainment on the calculation of the liquid holdup 


when valve V, is closed. To determine the mass 


flow of liquid drops in the cross-sectional area of 
the probe tip, the connection MC is blocked, V, 


is also discussed. 


EXPERIMENTAL TECHNIQUI fully opened and the water discharged through 


In order to separate the drop-laden vas core this valve collected and measured, To measure 
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the impact head, the 6in. length of the 
the end of the tip and the horizontal 
connection M¢ filled with fluid 
by purging with either air or water through RT, 
and V,. 


In arrangement 2. 


probe 
between 


i single 


must be 


whole stream 


Fig. 1, the 
emerging from the glass tube GT is confined in 
a polythene beaker B followed by 
bag PG: the liquid drains through WE to the 


feed tank and the air escapes through the hole 


a poly thene 


AE. The long probe PL can be used as before 
through hole H, or a shorter prob PS 4 in. long 
inserted through hole H. 

The manometer M is a U-tube 3 ft. long filled 
with bromoform (s.g. 2-9). The leg M, is open 
to atmosphere and M, is connected to the probe 
via a plastic tube PT kept in a horizontal plane. 
The bromoform in M, is kept under water which 
also fills tube PT ; the capacity of the horizontal 
portion of PT is such that it acts as a constant 
2-4). 


entrainment 


liquid reservoir (effective s. 


Measurements of and impact 
pressure were made on air and water flowing 
id. tube 


Continuous 


cocurrently upwards through a 0-427 in. 
for three different 


lengths of glass tube were used 10-25 in., 34-5 in. 


tube lengths. 


and 72in. long. Most of the experiments were 


Sampling devices. 


carried out with the probe in the centre of the 


core but some radial profiles were also explored, 


Limitations and accuracy of the experimental 


technique 


To obtain a truly quantitative sample of a 
discontinuous phase (water drops) flowing In a 
continuous phase (gas or vapour) it Is necessary 
to sample isokinetically ; only then are the stream- 
lines of the continuous phase directly in front 
of the probe undistorted by the presence of the 
probe. Sampling is isokinetic when the impact 
pressure registered by the probe is zero, since 
then both phases flow through the probe with the 
same velocity as they flow outside and upstream 
of it. 

In these experiments, since the path through 
the probs offers consice rably higher resistance 
than the 


kinetic. 


alternative one, sampling is not iso 


The smaller the probe resistance, the 


closer the sampling approaches isokinetic condi- 


tions ; this is why the shorter probe PS was also 
used, 

As the hydraulic resistance of the probe path 
inereases—by closing valve V, — progressively 
less air flows through it until at the limit, when 


the probe exit is sealed off, the gas streamlines 
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(Fig. 3). 


efliciency 


than 99 


dimensions 


under the 


conditions. 


Fie. 3. 


IS 


per 


hit the probe (see Fig. 2). 


hic. 2. Air purged probe 


related 
by 1 
Applying this plot to the 


and drop sizes of 100 and 50 4 


Gas streamlines 


Liquid drop paths 


The proportion of drops hitting the probe 
those crossing the same area in the 


the probe is termed the capture efliciency, ¢ 


suffer the greatest deviation. Under these extreme 
conditions the water drops, because of their high 


momentum, tend to cross the gas streamlines and 


to 
ol 
it 


to the licguad, gas and probe 


JUsSOURD and SHAPIRO 


per cent and 99-4 


conditions used for measuring 


[2 
present 
experimental set-up, for a gas flow of 30 lb hr 
the capture 
cent 
respectively when the probe exit is shut. Hence, 
the 
entrainment, capture efliciencies should be better 


cent, despite the non-isokinetic 


From a study of the volume—time 


Obedience 
4p. oF 


Solid probe (or liquid purged) 
Thin walled probe (air purged) 


D, 


probe diameter 
drop diameter 
collection efficiency 
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plots of the water collected the mass flow density 
of the entrained liquid (G,") is estimated to be 
accurate to 2 per cent. The impact pressure 


measurements are far less accurate. 


EXPERIMENTAL ResuLTS AND REPRODUCIBILITY 


When arrangement 1 was used, the vas cor 
was distinctly visible on discharge because of 
the drops it contained. It was quantitatively 
observed that the number of drops in the gas 
core increased with increasing gas or water rate 
and the size of the drops increased with in- 
creasing water rate but decreased with increasing 
gas rate ; the diameter of the drops was estimated 
to range between 50 and 1.000 nu. This arrange- 
ment was unsatisfactory for Impact pressure 
measurements since the amount of air leaving th 
system via FE rather than KE was unknown. 

The entrainment results for the top of the 
tube using three different set ups are shown in 
Fig. 4. Values of the mass flow density of the 
entrained liqnuid in the centre of the gas core. i.e. 
tha liquid sampling rate divided by the sampling 
area, are plotted against the mass flow rate of 
gas, with the liquid feed rate as a parameter. 
Despite the good control of the liquid and gas 
flow rates and of the liquid temperatur: mainly 
19-21 C and always between 18 and 22 °C — the 
reproducibility is on the whole 20 per cent, no 
particular set-up showing a consistent bias. The 
corresponding values of the impact pressures 
are shown in Fig. 5, where the re producibility IS 
of the same order. The log-linear plot is used 
because it can conveniently combine a thousand- 
fold range of variation of the independent variable 
with clarity of presentation. Results for the 
bottom and the middle of the tube are shown 
in Figs. 6 and 7 and radial profiles in Fig. 8, 
Fig. 8 provides the justification for taking th 
radial distribution of entrainment within the vas 
core to be uniform within 20 per cent of the axial 


value. 


VeLociry or THE ENTRAINED PuHase 
It is possible to calculate the velocity of the 
drops in the following way. 


The differential equation for the conservation 


of linear momentum flow is 
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\ir water tlow 


Entrainment rates at the tubs 
axis 7T2Zin, above inlet. Fic. 5. Air-water flow. Impact pressures at the tube 
axis 7T2in. above inlet 


lop, short probe, beaker lop, long probe, water purged, beaker 
lop, long probe, beaker Top, long probe, air purged, beaker 


© Top, short probe, air purged, beaker 


Pop, long probe _ knife edge 
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entrained drops G,, 


Jensity of 


‘ow 


Mass 


Aw rate W, b/w 
Fic. 6. Air-water flow. Series Wy, lb/hr. 
Entrainment rates at the tube B 
axis 32-5in. and 10-25in. above ( \ir—-water flow. Impact pressures at the tube 
inlet. D axis 


$2-5in. and 10-25in. above inlet. 
© Middle, short probe, beaker Kk 2: © Middle, short probe, air purged, beaker 


Bottom, long probe, beaker KF } Bottom, long probe, air purged, beaker 


x Bottom, short probe, beaker G Bottom, short probe, air purged, beaker 
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. My 

: 

be ous 
¥ tube . 
Fic. 8 \ir-water flow. Radial traverses at 7T2in. from inlet. 
(1 he (2) lb her (3) 220 Ib/hr (4) 1001b hr. 
Wy lb hr Mr 5-5 tb he 
Fores mass rate velocity Impact pressure of liquid drops (38 
or considering unit cross-sectional area in the where « is a constant depending on the behaviour 
dire t ol llow. ot the drops on reac hing the probe ; if the drops 
' follow the gas streamlines 2 1 2: if when they 
Pressure, P mass tlow densitv. G 
te ' hit the probe they are brought to rest and then 


For the continuous phase, the mass flow density 
net f the velocity ou) and (1 
tevrated between upstream conditions, 7 
real tive Vik lds 
impact pressure of gas pu* 


local velocity /mean bulk gas 


where A 


velocity 


gas 


Strictly, this impact pressure should be multi- 


plied by the factor R, R, R.. but this is 
virtually unity since R, R, (see below). 
For the discontinuous phase, (1) mtegrated 


over the same limits vields : 


escape along the probe face in a direction normal 
to their original path, G, 
and uf 


their velocity cl 


is ince yx ndent of the 


velocity they bounce backwards, 
inging sign but not magnitude, 
2. Suariro allows for the first two possibil 
third. 


mass flow of entrained 


ities but does not consider the 
Writing the ratio of the 


liqguid to the mass flow of vapour 


drop velocity u, 
local gas v« locity uy A 
(3) becomes 
impact pressure of liquid drops xtpSK —— (4) 
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T abl Calculations on interpolated entrainment and impact pressure val nirve of the tule during cocurrenl upward air wate flux 
Top Viddle Bottom 
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Thus the resultant impact pressure P, is given 


-{4K? +2) (5) 


The value of K can be obtained for fully developed 
flow in tubes from the universal velocity profile 
[3]. 

Calculations using entrainment and impact 
pressure values interpolated from Figs. 4 to 7 
and hold-up values from experiments previously 
reported [1] are shown in Table 1. The resulting 
values of a8 range mainly between 0-4 and 1.3 
but in extreme cases of very low values of v; 
values as high as 4-0 or more are obtained. Since 
it is the gas that entrains the liquid, S has an 
upper limit of 1-0 and therefore «S cannot be 
greater than 2-0. The abnormally high values of 
as may be partially accounted for by the esti- 
10° lb/ft hr? 


in the measurement of impact pressure). Since 


mated experimental error ( 15 


no concrete indication of the value of x is furnished 
by these calculations, its most probable value, 


a = 1-0, is adopted, 


SomE CHARACTERISTICS OF THE ENTRAINED PHASE 


A study of Table 1 indicates that the amount 
of liquid dispersed in the gas core is perhaps an 
exponential function of the distance from the 
inlet, first increasing quickly from a value of zero 
at inlet and then more slowly to a steady equili- 
brium value. However, if entrainment is primarily 
caused by the interaction of the gas core and the 
liquid interface due to the velocity difference 
between them, no true equilibrium is ever possible 
because the gas velocity continuously increases 
due to the unavoidable pressure gradient. 

The available data seem to indicate that for 
a given run, the process of entrainment starts by 


the detachment of a few drops which accelerate 


to a high speed and then, as the mass flow of 


liquid in the core increases, the mean velocity of 


the drops is reduced, 

From the radial traverses it seems that the 
mass flow rate of the drops and their velocity 1S 
reasonably constant near the axis of the gas 
core. As soon as the probe approaches the gas 
liquid interface, the ripples interfere with the 
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drops. Since the diameter of the probe is about 
20 per cent of the gas core diameter, the accuracy 
of the traverse is limited, but the use of a smaller 
probe would lead to inaccuracy on account of its 
diameter being comparable with that of the drops. 

No quantitative measurements were made on 
the size distribution of the drops, but it is possible 
to calculate the maximum diameter D,, of stable 
drops given the physical properties of the system 
and the relative velocity by the use of the relation 
originally proposed by Hinze. 


Critical Weber number — 


' 
(uy u “py 
constant. 

o 
Values of this constant have been proposed as 
20 and 18-2 for the higher limit and 11-2, 10 and 


For the air-water system, taking a value of 


8 for the lower limit 


(u, Up) = 50 x 10* ftyhr to represent condi- 
tions during the formation of the drop, the 
maximum drop size is 400 u. 


Errect OF ENTRAINMENT ON THE PREDICTED 
AND CaLcuLATED VaALuEs or Lieurp 
THICKNESS 


A previous paper [1] shows how the liquid film 
thickness, 
Y; pu*y a can be predicted. This prediction 


expressed in dimensionless form 
depends on the liquid flow in the film and so, 
when allowance is made for the proportion of 
liquid travelling in the vapour core, the predicted 
value of is reduced. 


The value of y; calculated from experimental 
results is essentially the product of the liquid 
hold-up R, and 4/7,. The value of R, measured 
is the sum of the liquid hold-up in the annulus, 


R,, and that in the vapour core, R 


e 


Since the liquid in the core travels faster than 
that in the annulus, its residence time in the test 
section is lower. If Uy is the fraction of liquid 


travelling in the vapour core, then 


R, = ad 
Ry Suy 


by 
VOL. 
l 9 
1960 
4 > > 
Le. Ry = Ry + R, 
|| 


i. 


rr) py (I 
RaW, 
OW, py 
Wy py 


R, ) 


alan al 


that 


to derive the 


It is 


from tha 


possible 
Aives [5). 
experiments Ss py 62.2 lb ft 
lp 13 ft® Ib. that the 


given as a percentage ol the total liq nd tlow. 


he Tht 


data of assum in hus 


and 


and entrainment ts 


146 0-051 O-O4S 0.000087 
1h 
146 


O.038 


0.00037 


0-000 0.00109 


It is seen that, even at quite high entrainment 
rates, 
to R, generally with little error, and R, l R,. 

Im the part the 
total 


R, and so Ry, can be taken as equal 


presence of entrainment 


pressure gradient is accounted for by the 


acecleration of liquid droplets and so the caleu 


lated value of will be reduced, resulting in 


Value ol y Is, when 


account ts taken of the presence of entramed liquid 


n lower calculated 


both the predicted and calculated values of y; 
are reduced. 

By a momentum balance, it has been shown [1] 
that in 


the absence of entrainment 


dP,, dP, aP 
iz d ad 


When entramed drople ts are present this equation 


is moditied to 


dP,, dP, 4P, 

d dz I: 
dW, 


if anv forces due to droplets leaving the vapour 


core at a higher velocity than ww are neglected, 


This can be written as 


a 


dP, dP, AP, 
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. dP, dP. dP, dP, 
ds dz d dz 


Lpy 


u, Gy ry) 


d 
(1 


G, 


dP. ry 
dz R, 


SG, 2 d | 
dz | 


Pa R, 


Evaluation of (7) and (8) from interpolated 
data as shown in Table 1 shows that there is con 
siderable improvement in the agreement between 
the predi ted and calculated values of Y, when 
entrainment ts allowed for (Fig. 9). In carrying 
out this computation the following assumptions 


were made: 

(1) The mass flow density of entrained liquid at 
the limits of the 
This 


anv cross-section within 


vas core ts equal to the axial value. 


has experimental support (Fig. 8). 


Top and bottom values of 7, and of 2S are 
used to evaluate the acceleration pressure 
gradient, assuming « to be constant along 
the tube. 


(3) When «8 1, S is taken as unity. 


Although the ambiguity in the corrections to 


U7) is large (of the order of 20 per cent) the 


ambiguity in the corrected values of y; is small 
that 3 The 
corresponding to both low liquid and low gas 


(Fig. 9) are left 


(generally less per cent). points 


rates unchanged because no 


entrainment values could be measured (mainly 
due to the resulting impact pressure being too 
the fluids through the sampling 


small to drive 


probe), 


CONCLUSIONS 


(1) The proportion of the total liquid flow 
the 


considerable at 


travelling as drops in core during 


annular flow can be very 


high gas rates; figures of 40 cont 


~ 


per 


uy R p, ae as 
putting rp and S 
9 
12 
1O¢ 
(2) 
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. 9 Comparison of theoretic al 


© original points 


entrainment (see Table 1, Series E) have 
heen recorded with indications that they 


further increase with increasing gas flow 


until virtually all the liquid flows as a 
dispersed phase 

The minimum obtained by the entrainment 
values in the case of the two highest licquic 
rates might be attributed to the transition 
from slug flow to mist flow ; nevertheless, 
the theory based on annular flow is obeyed. 
After the 
made, the experimental vaiues of y,; do 


the 


correction for entramment ts 


not deviate in general from values 


memento 


corrected for 
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and experimental values of ¥, 


entrainment. 


predicted from the theory discussed in 
detail elsewhere [1], to an extent greater 
than is attributable to experimental error [3]. 
Due to the absence of any guiding principles 
for the prediction of the extent, quality 
and distribution of entrainment, no correla- 
tion is attempted of the available data ; 
by inspection of Figs. 4 and 6 it is possible 
to say that for a given position in the tube 
the following is approximately true for an- 
nular flow 


conditions ; is proportional 


to W?* for constant W,, and proportional 
to W > for constant W, 


Y 
A 
/ 
4 
4 
9 a 
12 fy 
‘i A 
Y/ 
A 
| | 
expe 
(4) 
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NoraTion dynamic viscosity 
density 
Symbols 
D diameter 
¢ capture efficiency e. ratio of entrained liquid 


surface tension 


shear stress 


collected to that in the same cross-section in the Subscripts 


absence of the collecting device a acceleration 


superticial mass flow density A = liquid flowing in annulus 


actual mass flow density D entrained drop 

ratio of local gas velocity to mean bulk gas velocity f entrained liquid in vapour stream 
pressure i vapour liquid interface 

fraction of cross-section occupied by one phase I impact 

ratio of entrained drop velocity to the gas velocity L liquid 

velocity ! vapour or gas 


friction velocity at the wall : wall 
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Studies in optimization 


The optimum design of adiabatic reactors with several beds 


R. Arts 


Department of Chemical Engineering, University of Minnesota, Minneapolis 14, Minnesota 


(Received 1 December 1959) 


Abstract—The optimum design of an adiabatic reactor involves the choice of the size of each 


bed and the amount of interstage cooling between beds. These decisions can be made in an 


orderly way by applying the notion of dynamic programming, and in the case of a single reaction 


allow of a simple graphical presentation. The principal subject of this paper is case of a single 


reaction with costs proportional to the bed volume and the amount of heat removed. It is shown 


that this simple model suffices to solve a number of more general problems 


la 


de chaque lit et la valeur, a des stades intermédiaires, de la réfrigération entre les lits. Ces conclu 


Résume projet optimum dun réacteur adiabatique met en jeu le choix de la dimension 


sions peuvent étre tirees méthodiquement par lapplication de la programmation dynamique et 


dans le cas of une réaction unique permet une représentation graphique simple. Le sujet principal 


de ce mémoire traite le cas d'une réaction simple dont le coat est proportionnel au volume du 


lit et a la quantitée de chaleur ¢liminée. On voit que ce simple modéle suflit pour résoudre un 


certain nombre de problémes les plus généraux 


Zusammenfassung Die optimale Gestatlung eines adiabatischen Reaktors liuft hinaus auf 


die Wahl der Abmessungen jeder Kontakt-schicht und den Umfang der Kiihlung zwischen den 


Kontakt-schichten Diese Entscheidungen konnen durch Anwendung der dynamischen 


Programmierung getroffen werden. Im Falle einer einfachen Reaktion ist eine einfache graphische 


Darstellung moglich. Hauptgegenstand dieser Arbeit ist der Fall einer einfachen Reaktion, bei 


dem die Kosten proportional dem Kontaktvolumen und der abgefiihrten Reaktionswiirme sind. 


Es wird gezeigt, dass dieses einfache Modell geniigt, um eine Reihe allgemeiner Probleme zu 


lOsen. 


1. INrTrRopUCTION the composition and temperature variations are 


Tur ryre of reactor to be considered consists of related, for, if the reaction is exothermic. the 


one or more adiabatic beds of catalyst in which a heat of reaction is entirely used in heating up 


single reaction ts taking place. Before entering the reaction mixture. For a fixed reference 


the first bed and between each bed the reaction composition, from which the reaction mixture is 


mixture is heated or cooled to the most suitable obtained by letting the reaction proceed to a 


temperature, Such a reactor is commonly used, certain point, the composition at all times can 


as, for example, in the oxidation of sulphur be specified in terms of the concentration of one 


dioxide or the production of methanol. The reactant. Thus, for example, if the reference 


cooling between stages is needed to keep away composition of the reacting mixture were 


from equilibrium conditions and is accomplished O,: SO, : SO, = 34: 65:1 (per cent volume or 


either by an interchanger or by adding a cold mole fraction) then when the concentration of 


shot of reactants to the reaction mixture; in SO, is c, the stoicheiometry of the reaction gives 


this paper we shall be concerned only with concentrations of (0-345 0-33 c)/1-005 and 


interchanger cooling. (0-66 — 0-67 c) 1-005 for OQ, and SO, respectively. 


When a_ single reaction is considered the We shall use the concentration ¢ of the product 


design of such a reactor is susceptible to a simple — to define the composition of the mixture so that, 


graphical procedure. Since the beds are adiabatic if T is its temperature, the state of the reactants 
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RK. 


at any instant can be represented by a point in 
Let 


the rate of reaction be a function of composition 


a eT diagram such as is shown in Fig. 1. 
and temperature R (e, T), then if v is the linear 
velocity of flow and w& the distance along the bed, 


a mass balance gives 


de 
da 


R T) 


emperoture 
adiabatic paths in the concentration 


temperature plane 


Since there is no heat loss a heat balance over 


an element of the bed gives 
dT 


(=) 


H Rie T) 


is the specific heat and p the density 
H is the heat 


h and dividing 


where ¢ 
of the 


reaction. 


reaction mixture and of 


Writing Hc, p 


by (1) we have 


(2) 


dT 
di 


(3) 


h will generally be a function of c and T so that 
(3) 
has to be solved for the relation between ¢ and 


is a first order differential equation which 


For simplicity we will 
the 


equally valid for varying h) and the adiabatic 


T in any adiabatic bed. 


assume fh is constant (though analysis is 


paths are then straight lines with slope I h, as 


shown on Fig. 1. Thus if cy. 7’, is the state of the 


Ants 
reactants at the inlet to a bed (the point A in 
Fig. 1) the composition and temperature at any 
other point are related by 
Ty (4) 


It follows that at any point in the bed the rate 
of reaction Is a 
R (e, T, h (e 

(1) 


setting f r/v, 


only, 
T,). 


immediately, for, 


function of ¢ namely 
integrated 


the 


Co. Kquation 


can now be 


variables separate to give 


de 


(c) 


(5) 


where (f — t,) is the holding time of a bed that 
suflices to raise the concentration of the product 
from cy to c. The adiabatic paths can thus be 
graduated in the variable ¢ and the bed size 
required to get from a state A (cg, Ty) to another 
state B (Cy. T,) is the difference of the graduations 

The reactions for which this type of reactor ts 
used are generally limited by equilibrium and at a 
given temperature T the reaction cannot yield a 
greater concentration than ¢, (7°). If this function 
is plotted in the ¢ T diagram it gives a curve, 
R (ec, (T), T) = 9, 

hand side of this curve. 


C,, along which R and we 
are limited to the left 
Since the reaction velocity ts small for low tem- 
peratures and vanishes on the equilibrium curve, 
C,. there must be a poimt on each adiabatic line 
at 


which R is greatest. 


at which 


Let such a point be 


dR 


h 0, (6) 


the 


intersection 


and denote locus of these points by C 
Then the of this 


adiabatic line gives us a suitable 


curve with an 


origim for the 


integration in (5), and we can graduate the 


adiabatic lines with 


d 


7) 
Ale 


tic) 


Cu) 


Fig. 2 


reaction 


shows some ty pical curves for the reversible 
A= 
A and B. The curves C,, C,, and some of the 


~ B tirst order in the concentrations 


of 
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Fic. 2. Concentration-temperature plane for first 
order reversible reaction. 


contours of constant ¢ and R are shown; C,, 
by definition the curve ¢ 0, whilst CL is R = 0 
ort 0. 

The interstage cooling that is required to bring 
the reactants away from equilibrium can also be 
represented on the c-T diagram. Supposing the 
reaction mixture enters the first bed in a state 
corresponding to the point A (Fig. 1) and the 
bed is of such a length that it leaves it in state B. 
In the interchanger it is cooled down to a lower 
temperature but does not change in composition. 
This is represented by a horizontal line BC, 
where the abcissa of C is the temperature to which 
the reactants are cooled. In the next bed the 
reaction will follow the adiabatic path through 
C, reaching the state D, for example. Finally 
we may observe that as originally available to 
the reactor the reaction mixture will be of a 
certain composition but probably not at the best 
temperature ; the point O represents such a state. 
It is therefore necessary to heat, or cool, the 
reactants to the inlet temperature most suitable 
for the first bed. The decisions that have to be 
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taken in seeking the optimum policy for an 
original state O are thus connected with the 
points A, B, C ete. and we might expect that in a 
good design the line segments AB, CD. ete. 
would lie over the curve of maximum reaction 
velocity. To formulate the optimum problem 
precisely we must first construct a suitable profit 


function. 


Prorit FUNCTION AND 


MAXIMIZING 


Meruop or 


beds of an N-stage adiabatic reactor will 


be numbered backwards, N being the first bed 


the reactants enter and 1 the last. Primed 
quantities will denote exit values and unprimed 
inlet values. Thus c’, is the exit concentration 
from the last bed (and 7", the corresponding 
temperature) ; ¢, 1s the inlet concentration to the 
last bed and equals c’,, the exit concentration at 
the penultimate. We may consistently write the 
feed condition as T’y., and 
n N. Thus the total increase in 


If g is the tlow rate and « the value of a unit 
1 


quantity of the product, the rate of increase in 


value of the product stream is «q — 
This is the simplest expression for the rate at 
which money is made by the reaction process ; 
more sophisticated profit functions will be 
discussed in the last section. 

The cost associated with the process are those 
of providing and maintaining the catalyst beds 
and heat interchanges. We will again consider 
these as simple proportionalities, reserving 
generalizations until later. The cost associated 


with the bed will be taken to be proportional 
N 

to the total volume gq 2 t,. This includes the cost 

of the catalyst, of maintainance and of construc- 

tion amortized in the conventional way. The cost 

of heating or cooling will be taken to be pro- 

portional to the total heat added or removed, 

namely T,, 7 n+1l° 
1 
two unit costs the net profit from the reaction ts 


If 8 and y are the 


: 
: | #+0-0C 

\ 
08 \ 
\ 
| 
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If we write A= p yc, p/x then the 
total profit is qa Py, where 
Py 
— At, — |T, — T.-,|} (8) 
1 


The optimum problem may now be stated in 


precise terms Given the initial state 
and the constants \ and the bed sizes 
t,, to. ... ty and inlet temperatures 7, 7, Ty 
are to be chosen so that the greatest value of 


Is this 


somew hat 


attained. Although 
restricted 
later that it can be extended to cop 
complicated When thr 


the resulting maximum of P. will be 


secin a 


problem it will be shown 
vith more 
situations howe ol 


has been made. 
function 


only of Cy and 


(Cray Max P, 
It is thus necessary to make a choiwe of 2N 
variables and if these have to be found simul 
tancously the problem is formidable. Hlowever, by 


using the methods of dvnamn programming the 


probl be re clu ad probe This 


each of which requires the simultaneous choices 
of only two variables 

To show this we first observe that the \-stage 
process the profit function can be broken 
into two parts, the first stage (labelled NV) and an 
tN 1 )-stage process The exit onditions 
cy. Ty of the first stage are the feed conditions 
for the subsequent N 1)-stage process The 
optimum from the N-stage process with respect 
to its feed ey., will certan not he 
achieved unless the (N 1)-stage process 1s 
giving its best with respect to its feed ¢y, Ty. 
If we know the optimum (N 1)-stage policy 
we can get the optimum \-stage policy by varying 
the operating conditions fy T'» of the first stage 
only, always using the optimum (N 1 )-stage 
policy for the subsequent stages By selecting 
the ty, Ty for which Py is greatest we shall have 


\RIs 


246 


the optimum poli v for N-stages. In our notation 


this may be written 


fy = Max { Py + (ex (10) 


where the maximum is over variations of ¢, and 


Px CN Ty Ty 
and 
di > (11) 
ty Rle T Cnt 
ex 
Kquation (10 vives a simple method = of 
successively calculating fo ...fyi Jo of 
course, 
We are thus led to consider the optimum 
design of a one-stage reactor first afterwards 


making it the second stage of a two-stage reactor, 


and so The so-called principle ofl optimality, 


on which dynamic programming rests [1], has 


heen invoked to establish the basi equation (10), 


3. Tue Orrimem PeRPORMANCE OF A 
SINGLE Bep 

For a single bed we have 

A Max p, 
Max ic, T, (12) 

where the maximization ts by correct chores ol 
and 7, Suppose that T, has already been 
chosen, then simee ¢ ( we are already 


committed to the aciabatic path through 


r,) and have to choose to maximize Py 
llowever from (12) and (7). 

é “ge 

\ ae 7 \ 

dt, 


being taken along the adiabati 
the point 


T i/ must the retore lie on the curve R T) 


the derivative 


path If p, is to be maximum end 


\. which we denote by T° Moreover only the 


part of the curve R \ which lies to the right of 
(,, sof interest since otherwise dp, / would be 
during the 


Analytically the condition for the 


negative reaction, which is unco 
unprofitable. 


end point is 


R, R 


R T, =A 


i. 
VOL. 
19% 
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Consider now the choice of the best 7, that is For given c’, ¢c, the optimum policy for the 


the point A on the horizontal line ¢ C's ¢, choice of 7, and the resulting maximum profit f, 
Denote by 7, the integral is now seen to have a very simple structure. 
e, In Fig. 4 S, and S* are the points where the 
1, [ (1 A de (14) tangents at 7. and 7*, meet the axis J, 0. 

R Then if: 


where c’, is always chosen to satisfy the condition 
(13), so that J, is a function of 7, only. If 1 is at AC 
C in Fig. 3 the intersection of c¢ e, and I”, 
then Cy and I, 0. T, decreases (A 
moves to the left) I, will increase at first and then 


decrease for in the position 4’ B’, where B’ is 


the intersection of c,, and J”,, J, is certainly 


negative. It follows that J, as a function of 7, 


must be a curve such as is shown in the lower - 
part of Fig. 3. But 


(i) 7’, < S,. p, is always negative and it is 


unprofitable to react (i.e. the reactants are 


so cold that it is not worth heating them 


up at the cost of heating ,). 


Py Maximum when 


T - 
(i.e. the reactants should be 
heated to Tye before entering the reaction 


bed). hy I, (Ty) (T,* 


T T’, T*,. py is maximum when 


7, T 0 (8.0. the reactants should he 


introduced into the reaction bed immed- 


f, = I, (7",). 


diately). 


T,* <T, <S*, py is maximum when 
T, z*, (i.e. the reactants should be 
cooled to T°, before entering the reaction 


bed). f, =1,(T*,) — — T*)). 


Fic. 3. The construction of the integral J). 


the difference between the humped curve of J, 
just constructed and the V-shaped = curve (v) S* T’.. p, is always negative (.e. the 
pw |T, T’... The maximum difference will reactants are too hot to be worth cooling). 


occur when T, has a value of either 7,. or 7%, 
the point at which 4, d7, pe Or — mp. Moreover the maximum profit function 


Sy T's) I, (c's, T's), <T,° > (15) 
(C 9, T’,*) T*,), S* | 


ws 
7. 
12 
Fic. 4. flee, T’..) for a given c’,. 
1460 2 
Py I, 7 1 
(ii) 
WN — 
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K 


eurve (ec, ») bevelled off 
with its two tangents of slop rT 

In the held 
constant and the values of 7T*, and 7T,. deter 
mined. If this is done for 
rm, and as functions of «¢,. Thes 


plotted in the C, 


Is precisely the 


been 


last paragt iph Le 


Various we obtain 
may be 
T as a curves having 
and J’. rh 
policy for a single bed ts completely ven by the 
Fie. 5 If the feed 
fin Fig. 5 


two branches optimum 


eurves I, and J”, as in 
point lies to the left of (as 


Pia. 5. Typical optimum policies with on 


he at thre 


then react 


reactants 
and until R 
If the feed pont les between the two 


the poles is to 
(point B) 
(point ¢ 

branches of I’, the policy is just to react (the line 
DE If the feed point hes to the right of (=, the 
reactants are cooled to z° ‘ reacted 
the line FGH). Since the single bed will 


be re varded “as the two 4 it is 


and 
next 
second of 
the last situation that is important. The surface 


(c's T 


.) is shown in Fig. 6. 


The optimum yield surface. 


The 


show the structure of the policy but to caleulat 


fore onstruction has been to 


thy eurve 


Thus (14 


Pr, it ws best to proceed analyte ally 


is actually 


1 
\ | | =) de 


In both the mteerrals and J, the lint 
fs so that on 
alculated as 


k 


hes on the aciabati 


path J, cnh 


curve any 
function of the 
limit by from the 


Now 


lower integrating upper 


cording as T, 


the 
the se oar th 


and 


Then on anv adiabatic path we can find 


which J A 


and 


pomts at 


intersections of the curves / with 
the acdiabatr 

For the 
partial derivatives of 
stated here : the 
appendix. In 
that J 
for its end points and R, 


path. 

next step it is necessary to have the 
These 


proof of them ts given in the 


these 


will only bn 


formulae it ws understood 
is calculated using the optimum policy 


R (e,. T,). 


T, being 


‘ hose optimally. 


If C » T 2 lies outside the branches of the curve 
r, then 
2 (e's) \J, fa. 


poor ji. In the important case 


Tu PERPORMANCI 


Beps 


two adiabatic 


or SEVERAI 
feed 


In the case of beds with 


conditions c’,, T’, we have 
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so 
: (17) 
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We again assume 7, the inlet temperature to 
the first bed to have been fixed and ask at what 
concentration ¢ 4 the reaction in the first bed 
be stopped to give the greatest value of p, hi. 


T 2 Is line arly related to c’, by 


T, ha 2 9) 


so differentiating along the adiabatic path 


Setting this derivative equal to zero for the 


(20) and (22) we have 


and using 


maximum 


Thus 


where 


R, Rie, 


= 


R(e,. T,) R, (233) 
and the reaction in the first bed should continue 
to such a point that on cooling to the optimum 
inlet temperature for the second bed the reaction 
will continue at the same rate. This is a remark- 
ably simple eriterion and allows acurve I" ° to be 
drawn on which the exit conditions of the first 
bed must lie. To construct I’, we have only to 
take a point ¢,, 7, on the curve J*, and draw 
through it the line of constant reaction rate R 


until it again meets the horizontal line ¢ Cy. 


Fic. 7. Exit conditions for the first of two beds. 
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It is evident that rs must lie between Case and 
I’,. This is shown in Fig. 7. Analytically we 


find 7”, as a function of c’, by solving 


até )] R le (c’,)]. 


To tind the best value of the inlet temperature 


T, we proceed as before. The quantity 


(Co, 


as a function of 7, for fixed c, gives a humped 
distortion of the curve of Fie. 4 to which we can 
draw the tangents with slope uw and obtain 
f, precisely as before. This will yield two tem- 
peratures 7',, and 7%, to which the reactants 
should be heated or cooled according as is IS 
less than T*, or greater than 7”. 


I’, with two branches can be drawn in the ¢, T 


Avain a curve 


plane and the policy for a two-stage reactor 1s 


found immediately in the same wav as before. 


Fic. 8. Typical optimum policies with two beds, 


Fig. 8 shows the three different kinds of reaction 
path given by a reaction mixture initially too 
cold (ABCDE), at a 
(FGHR) and initially too hot (LMN PQ). This 


process will generate an optimum policy surface 


suitable temperature 


f, (c's, T's) lying wholly above the surface f, and 
therefore extends to some degree the limits 
S,, S*, without which it is unprofitable to operate. 

Analytically the condition for the optimum 7, 


simplifies in a remarkable way and if 7’, > T, 


| 
\ 
Ie T's 
12 
\ \ 
he, R, 
5 K 
\ 4 
‘\ 
| \ 
re 
\ 
\ 
\ 
\ 
Go 
re a 


. Aris 


(24) 
whilst if derivative 
\J, 
by integrating back from E 
until J, 
to give 


in the App nadix, 


T, < T’, the same Is 
2u. Thus I*, and I, can be determined 
along the adiabatic 
2A/A 


The detailed manipulations are given 


0 to give I*, and on until J 


Even more remarkable is the 
fact that the partial derivatives of /, (¢’,, T's) have 
the 


the 


same structure as before and in particular 


an adiabatic direction ts 


~(1 


where R, is the reaction rate at the optimum inlet 
T,) for the feed 3: T 

Because of this we can proceed to the three 
the 


derivative in 


dT"; 


conditions (« 2: 


bed optimum with equal case, Is curve on 


three beds 


the 


which exit conditions of the first of 


should lie is constructed from J, im same 
wav as I”, was constructed from see Fig. 7). 


When I”, 


the curves for which J, 


is known, and yy can be found as 
Oand 2u/Ar spective ly. 


5. GENERALIZATIONS 


Although the profit function ?’, may appear 
bn adapted 
Until 
now the parameters A and » have been regarded 
but it the 


to be too simpk to be realistic it cat 


to meet a variety of practical problems. 


as constants, is clear that 


optimum 


policy is a function of A and a. The maximum 


profit function is thus fy A, 


particular if mw is fixed the total bed volume 
N 


Vy =¢q Xt, will depend upon A. By repeating 


the calculation for various A the function Vy (A) 
be 
for A is the relative cost of bed volume and the 
tine 


may found ; it is a monotonic decreasing one 


higher this is the smaller will be optimum 
volume. If then a design were asked for to give 
a certain total volume V it would be the solution 
of the optimum problem outlined above with 
A = Ay so chosen that Vy (Ay) V. This might 
be a consideration of overriding importance as, 
for example, when the reactor is required to 
fit into some existing high pressure forging, which 
valuable to waste. Programmed for a 


is too 


modern computer such a problem as this would 


be very rapidly solved. If the interchangers 


had to be 


within a given volume it would be necessary to 


between stages also accommodated 
vary both A and » and to use data on the volume 


required by an interchanger removing a given 
amount of heat according to best design practice. 
In this case for each » there would be a value of 
A Al the 


and the final choice of ~ would be the one that 


satisfying volume requirement 


made fy Ty 13 A (ye), pe) the greatest. 


Again supposing A and yw are known we 


might ask given output 
q take 


account the cost of pumping the reactants at a 


for a design with a 


. 
Cy): Here we may into 


exit concentration 
We can design the last bed 
Taking 
back from the point ¢ c 


rate For any the 
= Cyn + 
first for its exit conditions must lie on J 
the 
I”, to I'*, we have the inlet conditions to the 
last bed Cy; E a Moving to the right along the 
line ¢ c, to the curve 
of the 


down an adiabati 


adiabati , on 


I". gives the exit condi- 
tions \vgain moving 
T, on I*, and so 


adiabatic line will cross the horizontal 


penultimate bed 
we find Cs, 
on, Some 
line through the feed ¢ Cp and this intersection 
the feed to the bed. 


Thus the number of beds to give the required 


will give conditions first 
output has been found and the volume of each 
cak ulated. 


the pressure drop through the reactor 


may be Using the accepted design 


practice 
may be found and hence the cost g (q) of pumping 
the reactants through at a rate q is found. Finally 


q is chosen to maximize fy (eys,, Ty.) — g (9). 


If the value of the product is not proportional 


but is a more complicated function, 


N 
q = 
1 
precisely the same technique may be applied as 
that 


(Cyey) 


above save we redefine the integral J as 


if v(ec) is differentiable. The curve I", is now 
(c) Ric. T) A 
restriction could be included in 1 


given by 1 specification 


(c), for in case 


e must lie between c, and c* we could set 
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v (ce) 0 outside this interval. More interestingly 
v (c) might take account of subs quent operations 
on the product stream such as the separation of 
the valuable product. If the optimum policy 
of operating the separating unit is studied so that 
v (c) is maximized with respect to this operation 
then the principle of optimality shows that we 
would find the optimum overall policy for both 
reactor and separating unit, 

If the cost of bed volume and heat interchange 
is not proportional to the volume or the heat 
removed the optimum problem ean still be solved 
by solving the sequence of equations (10), using 
the appropriate non-linear expressions in py. In 
this case, however, the simple structure of the 
policy will be lost and it will no longer be 
possible to draw the sequence of curves I) and 
r 


powers of a modern digital computer, 


a» The solution, however, is still within the 
It thus appears that the method of dynamic 
programming is well suited to a number of 
significant problems in the design of adiabatic 
reactors. It may be extended to the desivn of 
cold-shot reactors and to the case of more than 


one reaction, as will be shown in later publications. 


APPENDIX 


To calculate the partial derivatives of f, we consider 
the three regions defined by (15). Take first the region 


Ty < T*, in which 


A 
le 


f (c 


. 
c 


and Ry (c’; \ and cy 


In the expression on the right side of « , isa function of 


ce’, and T">- Thus differentiating the limits and under the 
integral 


1 Ah J; 


9 
= 


But since R, \ the first term vanishes and 


A) 
At MJ, (25) 
| R,| 


Similarly 


=AdJ (26) 
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d A 
h hi (27) 
| R,} 
In the region S, Ths 
| \ | 
l ck 
Here 
| R Me | R,) 
f 
"dk | de’, 


but \and AJ So that 


| R,| 


Clearly df,/dT", ja \J,, so that the same formula 


apply, it being understood that the path for J 


is chosen 


optimally \ similar calculation for the region 
T,* a S* vives the same formulae 


To establish the formula (24), let 7 3 T*, and 


R[c, T, + h(e 
— Ty) +f, T's) 


and observe that c’, and 7”, are both functions of T, 


However since T,, hic’, we have 


Then 


AJ 


(29) 


2 ; 
to get the second line we use (25), (26) and (28) and the 
R’, and 


final reduction depends on the fact that R, 2 
T’, the sign of the term 


AJ, In case T, 
u(T’, — T,) in FP, must be changed, this gives a p in 
the penultimate line of (29) and so )F,/)T, AJ, — 2p. 
To show how the derivative of f, simplifies so dramati- 
cally we will consider the case T’, > T*, (c’,). Here 


JOL. 
12 
Q6HU 
1 h Zz. (28) 
aT. 
MF, A | & 24 4 IC» oT, 
\ A IC's 
| | R,} oT, 
| R 1! | 
oT, dT", 
and 
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thi ‘ in this region this formula can again be extended. 


and the R in J, is R fe, 


T both depend on ¢’, and but sine Finally 


) | 
h té — = (34) 


is everywhere valid, so that the derivative along an 


widiabatic path of the optimum /, is related to KR, just as 
the derivative of f, was related to R,. 


NovraTION 


concentration of product 

specific heat of reactants 

concentration for maximum reaction range 
maximum profit from N beds 


at of reaction 
in which the first line is the derivative of 7,. Making use iT] 


of (31), (25) and (26) and collecting terms together nite. defined by equation (16) 
integral defined by equation (17) 
hid number of beds 
7 total profit from N beds 
profit from bed 
flow rate of reactants 
(c, T), reaction rate 
temperature 
holding time, 
total volume of N beds 
linear velocity of reactants 
value function of product 
distance from inlet of bed 


since the square brackets vanish. Of course unit value of product 


the second term is really Au. However a little more work unit cost of bed volume 


shows that in the form (32) the expression is valid when unit cost of heat interchange 


I’, < T*,, it being understood that J, is on an optimal x 
path. Clearly ¥ Cp i 
density of reactants 


x 


Sullix n denotes value in the n'" bed from the exit. 


A prime denotes conditions at exit of any bed. 
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The dynamic characteristics of a rotameter type flowmeter has been investigated 


experimentally by making frequency response measurements using floats of various sh ipes and 


materials in two different liquids. The experimental results are compared with calculations 


based on a simplified theory which is derived. 


Résume 


La caractéristique dynamique d'un debitmétre type rotamétre a été étudiée expéri- 


mentalement dans deux liquides différents en faisant des mesures de réponse de fréquence, et en 


utilisant des flotteurs de formes et matiéres variées. Les résultats « xpéerimentaux sont comparés 


Zusammenfassung 


aves des calculs bases sur une theort simplifies quia ete établie. 


Die dynamische Charakteristik yon Rotamessern wurde durch Frequenz- 


gangmessungen experimentell untersucht in zwei verschiedenen Fliissigkeiten unter Benutzung 


von Schwimmern verschiedener Gestalt und Materialien. Die experimentellen Ergebnisse werden 


INTRODUCTION 


PREVIOUSLY published studies of rotameters have 
been confined almost entirely to considerations of 
their operation under conditions of constant 
liquid flow-rate and the re lationship between the 
flow-rate and the height of the float in the tube. 

and [1] derived a correlation 
between the flow-rate and the rotameter scale 
reading, using a modified form of Bernoulli's 
equation. ‘This method involved the graphical 
determination of two functions, whereas the now 
generally accepted correlation, first proposed by 


ScHoEN BORN and CoLBURN [2], requires the graph 
ical evaluation of only the rotameter coetlicient 


K pp, in the expression . 


Px 
where @ — the volume flow-rate, 
AK» the rotameter coetlicient. 
A the miminum annular area between 
the float and tube wall, 
g the acceleration due to gravity, 
i the volume of the float, 
Py the density of the float, 
Pu the density of the metered liquid, 


mit Berechnungen verglichen, die auf einer eigens abgeleiteten vereinfachten Theorie beruhen. 
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This expression was derived by assuming that 
the flow between the float and tube was equivalent 
to that through an annular orifice. 

Using the same basic assumptions, Martin [3] 
derived a slightly moditied version of the above 
equation by substituting a function of the 
Reynolds number for K p. 

The work of Fiscuer et. al. [4] was concerned 
with the design of a rotameter float which would 
have a calibration independent of the viscosity 
of the metered liquid. 

[5] found that A, 
was very nearly constant for different values of a so 
This method 


was suggeste d as a means of correlating the cali- 


Using equation (1), Frrrsem 
called ** viscous influence number.” 


brations of rotameters metering different liquids. 
Virovec and RezApek, [6] by dimensional 
analysis, obtained an expression similar to that 
derived by Martin [3]. 
These authors all concerned themselves only 
with the ste ady state operation of the rotameter. 
MacMIL.an [7 


expression for the transient response of a rota- 


|, on the other hand, suggested an 


meter but this would be applicable only in cases 
of laminar flow through the annulus between the 


float and the tube wall. 
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Tue Dirrerentiat Equation or TH agreement with experimental results, so that it 

UnstTeapy OPERATION OF A ROTAMETER was necessary to devise a means of investigating 
The forces acting on a rotameter float, displaced the nature of the function /, 
from its equilibrium position, are An experimental study was made of the 
1) The upward forces due to the drag of the fluid terminal velocities of steel and duralumin floats 
i! the innulus the pre SSUTE diff renee ACTOSS and ste el and bronze spheres (up to lin. diameter) 
the float and the buoyancy of the liquid in cylindrical tubes through which liquids flowed 
(2) The downward loree duc to the we ight ol the varsous rates (eet Fable 1). It found that 


float the « xp rimental results could be expressed in the 
oat. 


form : 
The forees due to the drag and the pressure drop 
will be some function of the fluid velo ity in the 


annulus and the vertical velocity of the float. 
where volume flow-rate; minimum annular 
Thus a foree balance across the oat may be 
area between the float (or sph re) 
written : 

and the tube, 


M dth » velocity of float or sphe re in the tube, 


dt? “dt value of u when J 0. 


where M the mass of the float, 


| the linear velocity of the liquid in the 


Table l 


annulus. 


the height in the rotameter tube, 


‘loat Tul 


diameter diameter 


Shape of Float 


float material 


L 
time, (ft) (ft) iquid 


and the othe r symbols have the san sionificance wuneter teel Water 


Rotameter Steel Water 
as bye tore 
Rotameter Steel O-0625 00753 Water 
The force, F. a ting on a to torm Rotameter Steel Glycerol 


drag is usually expr ssed by an ¢ j tion of the Rotameter Steel 17 Glycerol 
form Rotameter Steel 0-0625 0753 Glycerol 
Rotameter Duralumin OGvO Water 
Rotameter Duralumin 0625 O717 Water 
Rotameter Duralumin O753 Water 
Assuming that the pressure drop across the Rotameter  Duralumin 0690 Glycerol 
Rotameter Duralumin 0625 Glycerol 
annulus may be « xpre ssed in the same form as Rotameter Duralumin 0625 : Glycerol 
pressure drop across an annular orifice, the force Rotameter Stee! O852 Water 
acting on the float, due to the pres difference Rotameter Ste O28 Ose Water 
across it. mav be written Rotameter St O28 0428 Water 
Rotameter Ste O418 Water 
Rotameter 0540 Water 
Rotameter Ste 0466 0502 Water 
Sphere St OS13 Water 
Thus it is to expected that the function Sphere Ste OS135 0425 Water 
f(V, dh dt) in equation (2) would be of the same Sphere Steel 0313 O418 Water 
1) Sphere Steel O469 O05” Water 


form as equations (3) and 
Sphere Steel 469 O592 Water 


By analogy with a body falling through an 
Sphere Steel 0835 Water 


infinite fluid, the obvious choice for a first attempt sphere Steel 0833 OO16 Water 
would be the substitution of the difference  Spher Steel 0625 0690 Water 
between the float velocity and the mean velocity Sphere Steel 625 O717 Water 
Sphere Bronze -0690 Water 


of the fluid in the annulus, for V in equation (2). 
Sphere Bronze 0-0625 OOT1T Water 


This substitution gave results which are not in 
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d 
J(u — u,) ft ‘see 


Ad 


In Fig. 1 (d;/ Ad) (u — ug) is plotted against V and 
it may be seen that for all combinations of float 
or sphere, tube and liquid B is constant and equal 
to 1-65. 

At the terminal velocity of a body in a liquid, 


a force balance gives 


u) (5) 


gi (py Pw) 


where ¢ is an unknown function of V and u. 
¢@ may be assumed to have a form similar to that 
of F in equations (3) and (4) with a function of 
V and u replacing V*, the function being such 
that it 


each combination of float, tube and liquid. This 


is constant, with varving V and u, for 


condition is satisfied by taking: 


u d 


6) 
Ad 


— Pw) = Kr | 


Comparing this expression with equation (1) 
which refers to the ope ration of a rotameter in 
turbulent flow Ay 1/K,? if u 0. 


K, has been shown [2] to be a function of the 
R 


However 


Reynold umber 


Fic. 2. @ Rotameter with water. O Rotameter with 


glycerol solution. Cylindrical tube with water. 

Cylindrical tube with glycerol solution. 
Reynolds Number Rep Vp,,Ad/u. Ky was 
calculated for the floats in cylindrical tubes and 
Fig. 2 


Re, defined as 


shows the relationship between and 
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\d 


The values lie on the lines 


A» against Re R 


THE Freevency Respownsi 


or A RoTaMerer 


work on bodies through a 
tril it 


the 


From the move 
clerive an 


the 


was 
kinetn 


rotameter tloat when displaced from 


evlindrical possible 


expression for forcees on 


ts equilib 


Substituting the ex) from 


rium position 


equation 6) into equation (2 thy ree balance 


he« omes 


where d, is the diamet tube mn 


the displaced or trat lloat 


i-h 


@ may also be expressed in terms of the equilib 
im the 
diameter of the 


rotameter tubs 
tribe 


float 
Let the 


rium position ot the 
using equation 


at this point bye then, 


whence equation Ss becomes 


and W. D 


ARMSTRONG 


M dh 
g dt? 


p...) (10) 


If the 
the tube then put 


angle between the taym red sides of 


d 2H] tan 


d, 2h tan 


Then if is small, H*® tan? and tan? @ may 


be neglected compared with 2d, tan @ and 
tan @. 


Substituting A, A,? 


Hd, 
YA tan? ./* 


. equation (10) becomes : 


dh 
dt 
i,* dh\? 
(ZA pb h tan dt 


Put H andh H h 


that A’ and Hl’ are small: then fil 


at long yx nods of oscillation. i.e. if 


, and assume 

P sin wt, 
w 1s small, 
neglected and (11) 


dh may be equation 


become s: 


d, 
pb tan a 
| P sin wt (12) 


From this equation the frequenes response of a 


rotameter has been calculated. 


EXPERIMENTAI 
The apparatus for determining the frequency 
response of a rotameter is shown diagrammatically 
if Fig. 3. The liquid was pumped from the reser- 


voir 4 up to the constant head tank B, whilst 


d, Ad A 1 1, Dx A rim (d, d, 

aor: gB (d, — 
dt 2g d,)* dt 

Vol 
< 
1O¢ 
ul K, 

: 

I? — h? 

2A 7A; d th 

KrArp,d lh\* 
{ 
2 Bt(d, — \ dt 

T 49 (py Pee! 

ony 


The 


the flow to the rotameter was controlled by means 
of valves V, and Vy. Valve 1) was an ordinary 
globe valve by means of which the basic flow 
rate through the apparatus could be adjusted ; 
and VV, was the valve by means of which a sinu 
soidal flow was superimposed on the steady flow 
through V,. The constant level overtlow device 
C, was adjustable to enable the total pressure- 
drop across the flow system to be controlled. 

Attached to the driving shaft of the eccentric. 
which was used to oscillate valve V,, was a trip 
mechanism. This was arranged to close an electric 
circuit, operating a recording pen whenever the 
valve was fully open, Le. when maximum flow 
was possible through the rotameter Twin 
recording drums each with 25 min. capacity were 
used, The second drum wis used to record the 
instant of making rotameter readings. 

The experiment was repeated for several 
periods of oscillation of from 5 min to 2 see for 
each level of tlow. Three similar series of runs 
were carried out using ftlow-rates which gave 
rotameter readings near the top, middle and 


bottom of the seale, while the amplitude was the 


same in each case, The freque ney re sponses of 


rotameters with floats made of stainless steel and 


duralumin (diameter 0-75 in.) were determined 


frequency response of rotameters 


using water, and using a 45 per cent solution of 


glycerol in water as the metered liquid. 


RESULTS 


From a series of pairs of readings of rotameter 
lloat position and time, it was possible to obtain 
the best sinusoidal wave to fit the points by means 
of a statistical analysis using the theory of least 
squares. From this calculation, the phase lag 
and normalized attenuation were calculated. 

In practice, however, this method was found 
to be very laborious, involving large numbers of 
arithmetical calculations and a more rapid 
graphical method was devised. By adjusting the 
time seale, it was possible to plot all the readings 
from one experin nt on a single set of axes and 
obtain one cycle of a displaced cosine curve. 
From this curve, the phase lag was read off 
directly and the normalized attenuation was 
calculated from the amplitude of the wave, and 
the equilibrium flows with the oscillating valve 
in its two extreme positions. The accuracy of 


this method was checked ith several cases by 


Period of oscillation, sec 
Fig. 4a. 
e wate 
| 
e reodr« 
15.45 | 
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e 
~ 
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comparison with results obtained from a statistical 
analysis. 

The frequency response curves for rotameters 
metering water and 45 per cent glycerol solution, 
using stainless steel and duralumin floats are 
plotted in Figs. 4 and 5 together with the theor- 
etically predicted values from the solution of 


equation (12). 


Disct SSION 


From Figs. 4 and 5 it may be seen that the 


response of the rotameter does not change 
float woter appreciably with the liquid tlow-rate ; and it 
was not possible, within the limits of experimental 
s "t error, to detect any consistent difference in the 
- on results at the three flow-rates used in each case. 
~ > Thus, while all the results are recorded for the 
X frequency response in those cases where water 
ae was the metered liquid, only the mean values have 
a been recorded for cases in which the glycerol 
x 50 solution was used. 
The curves of phase lag and normalized at- 
Fie. 4d, tenuation show that the simplified differential 
equation (12) provides a re asonably good approx- 
imation to the frequency response of a rotameter 
at longer periods of oscillation, i.e. at lower values 
of the angular velocity w. However, at short 
periods of oscillation the predicted response 
deviates increasingly from the ¢ xperimental values. 
This is due to the fact that, in deriving equation 
(12). the velocity of the rotameter float in the 
< Bell tube was assumed to be small compared with the 
SS liquid velocity in the annulus and the square of 
the float velocity was neglected. This approxima- 
tion becomes progre ssively less accurate as w 
increases. The deviation of values, calculated 
using this approximation, from the expe rimental 
values would also be expected to be greater, when 
the coeflicient of the (dh dt)? term was increased 
as, for example, by decreasing the density of the 
float, or by increasing the density and the vis- 
cosity of the metered liquid. 

Thus the predicted response of the rotameter 
with glycerol solution begins to deviate 
appreciably from the experimental results at 

0 longer periods of oscillation than with water. 


Comparing the corresponding experimental 
Fie. 5b. results for the rotameter response with water 
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and with the 45 per cent glycerol—-water solution © = mass of float 


with a viscosity nearly four times as great, it @ = volumetric flow rate 


» » » ol l 

was observed that the effect of viscosity was me Reynolds number 
linear velocity. 

Small letters 

the extent by which the viscosity alters the d 


relatively small. The r sponse is affected only to 


diameter 
rotameter coetlicient A», whereas, in the case of g = acceleration due to gravity 


laminar flow, the coeflicient of dh ‘dt is propor- h = transient height in rotameter tube 


tional to the viscosity [7] and the difference ¢ = tame 


u float velox ity 


“tween the two responses is more marked. 
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Capital letters Subscripts 
cross-sectional area R appertaining to a rotameter 
a dimensionless constant T = appertaining to a cylindrical tube 
a constant é equilibrium value 
force f = appertaining to a rotameter float 
height of float in rotameter tulbx mean value 


dimensionless constant 1 appertaining to the metered liquid 
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Porosity distributions in a fluidized bed 


P. J. Baxxer and P. M. Heertses 
(Received 5 October 1959; in revised form 18 November 1059) 


Abstract Some theoretical equations, regarding the porosity distribution in a fluidized bed 
of verv large diameter, have been developed, Measurements on this porosity distribution in the 
centre of a thuidized bed of glass beads and air agreed sulli iently well with the theory. In general 
three distinet zones in the centre of the fluidized bed could be distinguished. Near the sieve, a 
sieve effect zone was found, followed by a zone of constant porosity re whing up to the initial 
bed height at incipient fuidization ind above it a zone of increasing porosity for which a pro 
bility function has been found. The constants in this function have been correlated with some 
fluidization variables, such as fluidization velocity bed weight and particle size 

The bed expansion in all cases showed a linearity with the product of bed weight and thaictiza 
tion velocity 

Experiments with silica vel t polystyrene pearls instead of vlass beads had the same results 

From measurements throughout the whole bed and by connecting points of equal porosity 
so called isopores have been obtained, from which streaming patterns in the fluidized bed could 


be deduced. The influence of the bed weight on such isopores has been pres nted 


Résumé Quelques equations theoriques ont ctees développées pour la distribution de la poro 
sité dans un lit Muidisé de trés grand diametre 

Des measures de la distribution de la porosité dans le centre doun lit fluidise forme des spheres 
de verre et de lair étaient sufissamment conforme au theort 

Een général trois zones pourt nt étre distinguées dans le centre du lit. Prés du tans se trou- 
vait une zone, appellee la zone de refit du tamus Au dessus de cette zone, une zone de porosite 
constante, ce qui sétendait jusqu’éa la hauteur du lit a Pétat prethaicdise Au dessus de cette 
zone, la porosité du lit dans la troisiome zone, augmentee jusqu a lunité. Dans cette zone une 
fonction «ce probabilite a cte trouvce pours la porosite, Les constantes dans cette fonction ont 
étées correlées avec quelques variables de fluidisation, comme la vélovité, la masse du lit et le 
diamétre des particules (spheres 

L’expansion du lit changeait lineairement avec le produit de la masse du lit et la vélocité de 
fluidisation 

Des expériments avec du sil wel and polystyrene spheres donnaient les mémes résultats 
quaves les spl res dle verre 

Des mesures dans tout le lit et par l'emploi des points de porosité égaux, des isopores ont 
étés obtenus. La déduction des profiles de vélocité pour lair et pour les parti ules était possible 


L influence de la masse du lit sur les isopores a ete presentes 


Zusammenfassung—Uber die Porositatsverteilung in einem Fliessbett von sehr grossem 
Durchmesser wurden einige theoretische Gleichunger entwickelt Messungen tiber die Porositats 
verteilung im Zentrum eines Fliesbetts aus Glaskugeln und Luft stimmten befriedigend mit 
der Theorie tiberein. Im allgemeinen kann man drei getrennte Zonen im Zentrum eines Fliessbetts 
unterseheiden. In der Nahe des Bodens ist eine vom Boden beeinflusste Zone zu finden, auf clie 
dann eine Zone konstanter Porositit folgt, die bis zur anfiinhglichen Betthdhe bei beginnender 
Fluidisierung reicht. Oberhalb liegt eine Zone steigender Porositaét, fur welche man eine Wahr- 
scheinlichkeitsfunktion angeben kann. Die Konstanten dieser Funktion wurden mit einigen 
Fliessbettvariablen in Beziehung gebracht, der Fluidisiergeschwindigkeit, dem Bettgewicht 
und der Teilchengrdésse. 

Die Bettexpansion hiingt in allen Fillen lineal von dem Produkt aus Bettgewicht und 
Fluidisiergeschwindigkeit ab. 

Versuche mit Silicagel und Polystryrolpleren anstelle yon Glaskugeln ergaben dieselben 
Eryvebnisse 

Aus Messungen innerhalb des ganzen Bettes und durch Verbinden der Punkte  gleicher 
Porositat erheilt man die sogenannten Isoporen, aus denen das Strémungsbild des Fliessbetts 


agbeleitet werden kénnte. Der Einfluss des Bettgewichts auf diese Isoporen wurde ermittelt. 
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OBSERVATIONS on the behaviour of fluidized beds 
at varying fluidization conditions and the results 
a bed fl, 2] 
have lead to a theoretical concept regarding the 
This 


concept will be pre sented in this paper tovether 


of porosity measurements in such 


porosity distribution in a fluidized bed. 
with the measurements which serve as an illus- 


tration of the theory put forward. 


Tue Porostry Disrripurion A 


Fiuipizep Brep 

As can be observed in glass columms, a fluidized 
bed, under not too extreme fluidization conditions 
can be divided into at least two layers, a first 
relatively dense laver from the bed support up 
to the height of the bed at incipient thuidization 
and a second relatively dilute layer above this 
height. The 


consists of two phases, a continuous and a discon- 


fluidized bed in the first laver 
tinuous phase (see e.g. Toomey and Jounsroni 
[3}). The continuous phase is formed of a bed 
of a relatively high packing density, the discon- 
tinuous dilute phase is formed of gas pockets 
ascending through the continuous phase. The 
continuous phase has a density equal to that of 
the bed in the pre-fluidized state. The ascending 
gas pockets transport a certain amount of partich s 
to the second layer of the bed. Circulation of the 
particles takes place, the particles moving down- 
ward pass along the wall of the bed. 

Instead of the term gas bubbles for the discon- 
tinuous phase as used by Toomey and Jounsrone, 
Wicke and Heppen [4] prefer the term cavity, 
since in fact no phase boundary is present. The 
cavity is a gas pocket which grows at its top and 
shrinks at its bottom and thereby moves through 
the bed. The gas pocket will therefore have an 
ellipsoid form, with the long axis in the direction 
of flow. The pockets 
through the bed is very irregular as has been 
shown by Martitrson, Hour [5}. 
The same authors have found that the increase 


movement of the gas 
and 
of the dimension of the gas pocket changes 
linearly with the height in the pre-fluidized bed. 

Interesting work has been done by Yasui and 


Jonanson [6] on the determination of the size 


characteristics of such gas pockets or cavities. 


Porosity distributions in a thuidized bed 


It was found by these authors that the length of 
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the pocket in the direction of flow increased with 
particle size, with the distance above the bed 
support and with the gas velocity. This growth 
with increasing height above the bed support ts 
mainly caused by coalescence of vas pockets, 
accom pank d by a decr ase in gas pocket frequency . 

In the fluidized bed of the first layer now two 
gas streams may be distinguished, a gas stream 
consisting of rising gas pocke ts and a gas stream 
through the continuous phase. 

From the above, it follows that the point 
porosity in a certain fixed place in the bed is 
never constant; at a certain time the porosity 
is equal to that of the pre-fluidized bed and shows 
the lowest value possible, some time later a gas 
pocket passes the point considered ; the porosity 
maximum value l. 


will then approach its 


Intermediate porosities will occur at other 


moments. These changing porosities are not easy 
to handle. 


an average 


therefore in the work to be described 
value has always been taken by 
integration over the time. Such a time average 
porosity appeared to be constant under constant 
fluidization conditions. This measured average 
porosity ts determined by the frequency of the 
vas pockets at the point chosen, with the restric- 
tion that the frequency and or sizes of the gas 
pockets is not such that overlapping occurs. 

If the time average porosity is assumed to be 
constant over the diameter of the fluidized bed 
an assumption which is close to reality if the 
influence of the wall can be neglected, for fluidiza- 
tion in which the gas velocity is uniform over the 
which 


chamber the arises as to 


distribution can be expected in the 


question 
porosity 
direction of flow. 
If the superficial gas velocity is V, a gas 
volume G, A V* passes through the bed per 
unit of time. This gas volume as said consists of 
two parts, the amount of gas necessary to main- 
tain the continuous phase, G,,, and the amount 
forming the gas pockets (G,). The amount 
G, is equal to the gas volume necessary for 
minimum fluidization and equals V,,, 4’ in which 
A’ is the cross-sectional area of the tube, not 


taken up by the gas pockets. This picture shows as 


*See nomenclature 


4 
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a consequence that an amount of particles 
vas pockets 
und of the 


is present 


corresponding to the volume of the 
present in the bed up to the height L, 
concentration in the pre-fluidized state 
in the second layer of the bed above the height 

The density p of the continuous phase in the 


first layer is determined by 


Pe (1 my 
if p, is neglected as compared with 


porosity at minimum fluidization veloc 


TRAWINSKI [7] has derived the sam type of 


equation for the viscosity : 


"Ne (1 (2) 


in which ec’ is a constant with the dimension of a 
viscosity. 

According to Marurson et al. 15) the velocity of 
the gas pocket will increase upon its rise through 


It will 


Ve low ity 


the bed, due to the growth of the pocket. 
be assumed that the mean gas pocket 


py) and 7 


V, is proportional to d,*, (p, 
The mean residence time ¢, of a gas pocket in 


the pre-fluidized bed is 
| 
t, (3) 
The total volume of cas pockets in tha first layer 
of the fluidized bed (from the bed support to the 
height L,,,) is then given by 


Vol, t, 


Since G = G, + G,y, 
V A’ and A A’ 


therefore VA G, 
G,/V,, it follows that: 


mf * 


A(V Lee 
| ms) (4) 


\ ol, 


As said an amount of particles corresponding to 
this volume must therefore be present above the 
pre-fluidized bed height L,,,. 

The porosity of a mixture of solids and gas is 


detined by 
Volume of gas 
Volume of gas and solids 


M 


therefore : L.A.p, 


in which M 
height. Thus for the pre-fluidized bed : 


M 
Lin A bs 


mass of the thuidized bed, L bed 


The volume of the solid particles in the bed above 


the height Ling: Vol, , according to equation (4) is: 


mi 


—V—) 


» mf 


(i 


Vol, 


The volume of the particles in the bed from bed 
support to Ling: Vol, is: 


M 
Vol, tmp (1 € 


mf 


The average y of the bed from bed 


support to L,, will be given by : 


porosity « 


or also by substituting (6) into (7): 


mf ( | 


“d Emp) 


mf} 
The variables in this equation are V and V,, 
because ¢,, and V,,, are defined by the system and 
conditions chosen. 

For a constant bed weight and particle size, 


equation (8) simplifies to: 


Emp + Vs) (9) 


For a constant fluidization velocity and particle 
(see above) and d, :: M 


[5] and Yasul 


size assuming V, :: d,? 


(according to and 


JOHANSON 


€d € mf weer V (10) 


mf 


The constant porosity «, from the bed support 


to the prefluidized bed height L,,, only occurs for 


mf 
points in the bed where restoration of the original 
situation after passing of a gas pocket can take 
place. Irregular behaviour can therefore be 
expected near the wall and near the bed support. 
Near the wall, because of a decreasing gas velocity, 
a decrease, near the bed support an increase in 
porosity are to be expected and have indeed been 


found. The question as to which height in the 
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= 2 (6) q 
€ mf Ith the 
a 
VOL 
“d 1 4 V €mp) (4) 
“mf ** Ps mf? 
Limi : 
G 
qr 
(5) 


bed the influence of the bed support will be felt, 
cannot be answered in advance. 

The particles in the bed above the pre-fluidized 
bed height L,,, have been thrown up by action of 
the gas pockets. The number of particles thrown 
up is a function of the number of particles around 
the gas pocket and of the size of the gas pocket. 
The height to which a particle will be thrown up 
depends on three factors : 

(1) The velocity of the gas pocket at the 

height L,,,. 

(2 The geometrical location of the particles 
on the gas pocket; a particle on top of 
the gas pocket will receive a larger impulse 
than a particle on the side of a gas pocket. 

(3) The number and type of collisions between 
the particles in the region above Ling: 
relative for example for co-current and 
counter-current moving particles. 

The distribution of the number of particles and 

therefore the relation between the density (1 €) 
in this region and the height above L,,, will be 


given by a probability function : 


(1 — C, exp | C,(h — Lyy)**) (11) 


For h = L,y, « should be equal to «4, therefore : 


(1 — «) = (1 — eg) exp [— C,(h — L,,)%] 


In the theoretical considerations as given above 
no influence of the wall of the bed has been 
incorporated, In general this will be permitted 
for columns with large bed diameters. In the 
experiments to be described a rather small bed 
has been used, influence of the wall will therefore 
occur. A number of authors [3, 4, 5, 8, 9, 10, 11, 
12, 13] have observed the downward stream of 
particles along the wall of the fluidized bed, at 
the bed support these particles are taken by the 
gas stream and carried upwards in the bed. 
Porosity differences between wall and centre will 
therefore undoubtedly exist. Fluidization vari- 
ables such as fluidization velocity, mass of the 
bed and particle size will greatly influence the 
porosity distribution in the dense phase (under 
L,,) and in the dilute phase (above L,,,). 

At high velocities no sharp distinction between 


dilute and dense phase will exist. 


Porosity distributions in a fluidized bed 


EXPERIMENTAL RESULTS 
The experiments on the porosity distribution 
in a fluidized bed have been carried out with 
glass beads of various sizes and with air as the 
fluidizing medium in a metal tube of 9 em dia- 
meter ; the bed support consisted of an electroly- 
tically made nickel gause with circular openings 


of approximately 80 « diameter. To approach a 
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Fic. 1. Porosity as a function of height in the 
fluidized bed at various velocities (centre of the bed). 
M 1.000 © 22-7 cm/sec 


d 175 210 


cm /sec 


O 59-6 em /sec 
78-8 cm/sec 


@ sec 


incipient fluidization (8-1 em /sec) 
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and 


s aA 


good, that is uniform, velocity protil the mcoming 
bed 
volume of 
rivid nickel 
place 


air passed a layer of Raschig rings under the 


small condenser with a 


support, \ 


approximately connected to 


tube could he placed in the bed a ny 


Full details of the 
elsewhere l 14 


wanted. apparat have been 


The t\ of the 


piven 


condenser was measured = with capacitive 


displacement meter, the resulting variable signal 
bration line 


bed 


measure ad. 


was integrated and by means of a « 
2 the mean porosity of the volume of the 


contained in the condenser could b« 


The 


tion of an alien body 


method chosen necessitated the introduc- 
into the fluidized bed. It ts 
quite conceivable that this would create distur 
pattern. It 


integrating the (« h) 


bances in the general flow has, 


however, been found by 
curves that the mean porosity value obtained in 
this differed in the than 
(maximum 8) per cent from the calculated porosity 
based on the height of the bed. 
The will be 


in graphical form in the next Figures. 


way mean less 


results obtained given manly 
For a bed weight of 1000 ¢ porosity measure- 
centre of the 


with 


ments in the column have been 


carried out fluidization velocities, varying 


from the incipient fluidization velocity (8-1 cm/sec 


of thundization velocity 


Ihe 


lines 


175-210 to OS-4em see, 
The 
obtained indicate the existence 

bed of thre« the 
be called the sieve effect 


for particle siz 
form of the 
in the centre of the 


results are given in Fig. 1. 


first which will 
zone (1) the porosity 
diminishes of heiwht in the bed. 
This is probably due to the lifting effect of the 


Zone, 


with imerease 


air stream entering the bed, as has been said 


Zone Il 
Between zones IT and IIL a transi 
found In zone ILL the 
For 


Il and in zone I 


before. is characterized by a constant 


porosity 
thon 


Is porosity 


increases vradually to unity increasing 


velocities the (h «) line in zon 
From Fig. 1 it can be 


Il into zone III 


approximately occurs at the bed height L_,. The 
PI anf 


moves to the « l 


seen that the 


transition of Zone 
(« A) line for Incipient fluidization for zone I 
differs from the others in this sense that « decreases 
near the sieve. This is probably caused by a 
compression of the pre-fluidized bed near the sieve 
by the condenser. 

The porosity «, has been plotted in Fig. 2 asa 
function of the fluidization velocity. 

In agreement with the theoretical considerations 
re lationship exists. At 


(equation 9) a linear 


incipient fluidization velocity the porosity of the 


bed should be equal to «,,, the porosity of loosest 
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packing possible. For a monodisperse system of 
glass beads, this should be the cubical packing 20 
0-476, which is remarkably tel 


with a porosity €m/ 


close to the value of 0-48 as read from Fig. 2. 


It has to be observed, however. that the system is ‘ 


not monodisperse, which creates a difficulty as to 


as to thy 


Porosity distributions in a fluidized bed 


aus a function « 


/ clearly, the type of sieve used and an 


j type of packing. 
| j Zone | 
L£ 

| 

tf effect caused 

al | | | / ol the Aone, 


yA / \ discussion ol 


fluidization 


the actual minimum fluidization velocity as well 
vibrations 


which are almost always present will disturb the 
discussion, 


in packing have 


large influence on the form of the lines and height 


later 


Fic. 3. The influence of the bed mass on the porosity 


distribution in the centre of a fluidized bed 


dy = 175 


cm see 


bed mass 


entrance 
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| 


Fic. 6 «q #8 a function of particle size. 


In another series of experiments the influence 
of the bed weight on the porosity distribution has 
been investigated. Results for bed weights vary- 
ing from 200 and 2000 ¢ at a constant fluidization 
velocity and again in the centre of the bed have 
been presented in Fig. 3. It can be concluded that 
no essential change in form and type of the lines 
occur at variations of the bed weight. The height 
of the zone II increases with increasing bed 
weight. For small bed weights (200 and 250 @) 
zone II disappears, only a transition point between 
zone I and III remains. The transition zone 
of II into III] however, increases in importance. 

From Fig. 4it can be seen that the bed expansion 
AL = L — L,,y, follows a linear relationship with 
the product of bed weight times fluidization 
velocity. 

A rather small influence of the bed weight on 
the porosity «, (or « of the transition point) 
exists, as follows from Fig. 5. If the assumptions 
giving equation (10) are valid the value of «, for 
large bed weights should approach «,,, 0-476. 
This, however, is not the case, e.g. due to slugging 
effects at larger bed weights the curve shows a 
minimum, whilst the value of « 0-476 is not 
reached. 

A third variable, the particle size d,, has also 


been incorporated in the experiments. For a — 
Fic, 7. Isopore-diagram 


1,000 g 59-6 sec 
= 175 — 2104 


constant bed weight (500g) experiments with M 


particles of different size have been carried out dy 
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M:500g 


M: 250g 


at various fluidization velocities. These are 
presented in Fig. 6, where «, is plotted against d,,. 
The term d, (max) from Fig. 6 refers to the 
maximum particle size, for which incipient flui- 
dization exists at the fluidization velocity used. 

So far, only measurements in the centre part 
of the fluidization column have been presented. 
The porosities out of the centre of the column are 
different from those in the centre. A_ radial 
uniform porosity distribution in the bed has been 
found. Therefore only the porosity distribution 
in a plane through the axis of the bed has to be 
measured, The results will be presented in the 
form of isopore diagrams, in which, in a plane 
though the axis, lines are shown connecting points 
of equal porosities. Such lines have been named 
isopores. An example of such an isopore diagram 
is presented in Fig. 7. From the isopore diagram 
streaming patterns in a fluidized bed can be 
deduced (see (1]). The presence of zones in the 
bed beneath L,,, near the wall of relatively high 
density indicates the existence of a whirl in the 
dense part of the bed. Variation of the fluidiza- 
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dp=175- 210 
V = 407 cm/sec 
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The influence of the bed mass on the isopore diagrams 


tion velocity does not change the general form of 
the isopore diagram; the influence of the bed 
weight, however, is notable. This can be seen in 
Fig. 8, where for a constant fluidization velocity 
the schematized isopores are given for bed weights 
varying from 250 to 1,000 g. It can be concluded 
that for large bed weights the porosity at the wall 
decreases. For small bed weights two whirl-like 
motions occur instead of one at larger bed 
weights. 

Besides with glass pearls some experiments have 
been carried out with silica gel of 5-5 per cent 
water content. Fig. 9 shows some measurements 
for various bed weights ; the type of diagram is 
the same as has been found for glass spheres. 

Finally the results as obtained for zone III 
(see Fig. 1) have been analy sed. The values of the 
constants C, and C, from equation (12) have been 
determined by trial and error in such a way that 
for each line a C, value has been established for 
which the value of C, was constant for each point 
of the line. The values of C, (1 — e,) derived 
from Fig. 2 and of C, and C, as determined are 
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Silicage!l 5.5 */e 
V2 75 cm/sec 

dp 7125-200 
@:-M: SOgram 

x M-:-100 gram 

o =150 gram 

 M=: 200 gram 


Fic. 9. Experiments with silica. 


viven in Table 1 for the 6 different velocities of 


fluidization for a bed of 1,000 g with particles of | besides «, apparently also C, ts a function of the 
175-210 


iluidization velocity V. In the region concerned 


It can be noted that C, is constant and that 


Table 1 Table 2 


19 
4 
A 
“ 200 
| 4 
11 
| 10}- ; 
ws (em) 9 VOL. 
xd 
2 
01 02 03 O04 OS O06 08 10 
M (g) 1 | Cs Cs, 
fem sec) ( 1 €d ( ( 3 j 
8-1 0-525 250 0-37 | 1-4 13 
22:7 0-46 0-59 13 500 0-42 | 13 
59-6 O15 13 1000 Ow O24 1-3 
78-8 0-30 O11 13 1500 0-37 O14 1-3 
O25 0-08 13 2000 | O36 1-3 
| 
im 
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from table I 
caiculated from 


(Vv Vnf)=78 


10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 


me  V (cm/sec ) 


Fic. 10. C, as a function of fluidization velocity. 


e@ Co from table I 
x CQ calculated from 


(C2-0.01)(M-110) = 220 
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C, as a function of bed mass. 
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and for a bed weight of 1,000 ¢ of particles of The relation between C, and d, can be re- 
175-210 yw the relation between C, and V can be presented by the equation (Cy - 0-8) 1500 
Cy,=78 (d,—d,__), which in d is the size of the 


min Pmin 
(see Fig. 10). particle just not pneumatically transported under 
In the same way the influence of the bed weight — the conditions used (see Fig. 12). 


expressed by the equation ( V | 


mf 


has been analysed, as presented for the same 

particles and a velocity of 40-7 cm ‘sec in Fig. 3. CONCLUSIONS 

Table 2 gives the results. From the experimental results as given it can 
Again Cs is a constant and equal to the value be concluded that the porosity of a tluidized bed 


found before, C, in analogy is a function of M. is not constant and the porosity distribution is not 


The relation between C, and M in approximation simple. However, a remarkably regular behaviour 


is given by : has been found for the time average values of the 
(C, — 0-01) (M — 110) = 202 (see Fig. 11). porosity. The equations as developed are covered 
7 rather nicely by the experiments. Rather large 

The same procedure can be applied to the variations in fluidization velocity and bed weight 
influence of the particle size d,. The results are do not affect the general pattern of the porosity 
given in Table 3. Again C, is a constant and (, distribution to a great extent. Even the dis- 
depends on d,,. appearance of the zone of constant porosity at 
small bed weights does not influence the validity 


Table 3 


of the equations developed, 


It can be expected that an increase in bed 
diameter will diminish the influence of the wall 
on the phenomena as presented, The porosity 
distributions as given in the centre of the bed will 


in this case be valid for practically the whole bed. 


NOTATION 
cross-sectional area of the fluidized bed = em* 
part of cross-sectional area of the 


fluidized bed cm? 


constant yoocm sec 


constant 
onstant 
onstant 
(1 €q) no dimension 
constant 
constant no dimension 


characteristic diameter of gas pocket em 


diameter of particle em 
volume of gas pockets per unit time em? ‘se 
volume of gas for minimum thuidization 

per unit time em? 
total volume of gas passing the 

fluidized bed, per unit time 

height in the fluidization column 

total height of the fluidized bed 

total height of the fluidized bed at 


incipient fluidization 
bed weight 


it 
“Oo 200 366 40c mean residence time of gas pocket 
in the bed 


Fa. 12. Cy as a function of particle siz fluidization velocity 
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Porosity 


mean velocity of gas pocket in the 
fluidized bed 
minimum fluidization velocity 


cm sec 
cm sec 
total volume of gas pockets in the 
fluidized bed 


volume of solid particles in the bed above 


em} 


the height Ling em? 


volume of solid particles in the bed from 


J 
em 


bed support to Leng 


porosity of the fluidized bed no dimension 


distributions in a fluidized bed 


€a porosity of the fluidized bed in zone II 


no dimension 


“mf 
fluidization 
density of the continuous phase 
density of the fluidizing gas 
density of the particles 
viscosity of the continuous phase 
viscosity of the fluidizing gas 


bed expansion 
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A bstract— WILKINSON and Arms Tron [1] have derived an expression which predicts the response 


of plate compositions in a distillation column to step changes in feed composition. This theory is 


valid up to moderate \ salues of tine but cannot be used to predict the changes in composition 
when the column approaches the new state of equilibrium. Following the same approac h it is 
possible to derive considerably simpler expressions which are valid for the initial part of the 


response. These solutions can then be extended so as to give the way in which the column ap- 


proaches equilibrium rhe theory is compared with experimental results obtained from a pilot 
distillation column 6 


Résumé— WiLkrNsown ct Anwerrone [1] ont établi une expression qui prévoit la réponse des 
compositions du plateau dans une colonne a distiller pour une variation pat stades dans la com- 
positions du plateau dans une colonne & distiller pour une variation par stades dans la composition 
de l'alimentation. Cette théorie est valable jusqu’aux valeurs modeérées du temps mats ne peut 


étre utilisée pour prévoir les variations dans la composition quand la colonne tend vers le nouvel 


état déquilibre. En attaquant k probléme de la méme fagon il est possible det iblir des expres- 
sions excessivement simples qui sont \ ilables pour la partie initiale ce la réponse Ces solutions 


peuvent étre étendues de fagon a indiquer la maniére dont la colonne tend vers l'équilibre. La 


théorie est comparée avec les résultats expérimentaux obtenus a partir d'une colonne de distillation 
pilote 


Zusammenfassung——WiLkrsson und AnmsTrone [1] haben einen Ausdruck abgeleitet, der die 
Antwortfunktion der Flissigkeitszusammensetzung auf dem Boedn einer Destillationskolonne 
ingibt fiir eine sprunghafte Aenderung der Zulaufzusammensetzung. Diese Theorie ist bis zu 
nicht zu grossen Werten der Zeit giltig, sie kann aber nicht benutzt werden, um die Aenderung 
der Zusammensetzung vorauszusagen, wenn die Kolonne sich dem neuen Gletk hyewichtszustand 
niihert. Nach derselben Methode kann man einfachere Ausdriicke ableiten,die fiir den Anfangsteil 
der Antwort giiltig sind. Diese Lésungen kOnnen entsprechend dem Weg der Anniherung der 


Kolonne an den Gleichgewichtszustand erweitert werden. Die Theorie wurde mit experimentellen 


Ergebnissen aus einer Technikums Destillationskolonne verglichen. 


INTRODUCTION section of the column shown in Fig. 1, following 


Twis paper is an extension of the work of & Step change in feed composition Arp gives 


Wi_Krxson and Armstrrone [1] for predicting the 
transient behaviour of a plate typ distillation Lip Oly 4 br,) + V (dye, oY,,) H 
column to changes in feed composition. The same 
theoretical model is considered and the theoretical where H is the molar plate hold-up, and dz, is 
basis of some of the steps taken in the present the displacement of the composition of the liquid 
work may be found in their paper on plate n trom equilibrium conditions. 
From the definition of the Murphree plate 
TREATMEN' efficiency based on liquid phase compositions 


A material balance over plate » in the enriching E = (z, ry) (@,* r,,) and assuming that 
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D X A exp nu, B exp nu, 0 n<N (4) 


n 


where 


(Sp 1) 
S; 


2p 


(seca) 


Similarly by considering a material balance 


over plate (N m) in the stripping section of 


the column it may be shown that 


( exp[(N Dexp| N mM) Vs 


lj 
m M (6) 


and ay 1s the slope of the equilibrium line for the 


Fic. 1. The column. 
stripping section. 


The boundary conditions requred to determine 
the functions A, B, C and D are: 


the equilibrium data for the enriching section may 


be represented by y Xpr,,* jp it may be 


shown that (7) Condenser 


Lp = Neglecting the condenser hold up the boundary 
oY, . [ (1 E) dr, 2) 


condition at the top of the column is 


E + E) 


Substituting (2) into (1) and transforming the 


equation by the Laplace transform gives 


1 (1 + Sp + p) X,, SeXnvi o 8) Taking the first two terms of the Taylor expansion 


where for X, gives 


a! x. (8) 


dn LE 


and the non-dimensional time 


apV(L—E) + | | 


Substituting (4) into (8) 


{ l 
2 Lp } 9) 


Following WiLkinson and Armsrrone [1] and 
reducing the finite difference equation (3) to a From (5) and (9) it may be seen that A and B 
differential equation by Taylor's theorem, the are of the same order of magnitude, for the special 
case Sp Lp 1, A 


general solution of XY, is 


| 
a 
| ly 
5 
Nor he re 
VOL. 4 J 
» 
| | ‘in —1 
/| (Ss Ys ~ 
+1 rg Se (1 + Sz) 
“Ss , 
EL, + (1 — EB) 
Ret er 
= 
273 
. 


Reboile r 


(ti) 


Similarly by taking material balance over the 
reboiler, the boundary condition at the base of 


the column may be shown to be 


| AX 
d(N + 
Ky 
pw a, S 
(10) 


Zs Sy 


where w is the ratio of the reboiler hold-up to 
the plate hold-up. Substituting (6) into (10) 


Dexp(N + M + 1)(v, — ®) 
agSp) E[ (as 1)S_/ag]+2, (11) 
(pu } E| (as l Ls Ue 
From the definition of v, and v, in (7) it may be 
seen that 
by | D exp aN + M 1) 


(iii) Plate N 
From a material balance over plate N in the 
unsteady state 
Ly Viyy., — Ledry 
\ 


putting n = N in 1 and subtracting this equation 


from (13) gives 


OYUN +1 BYN (14) 


The superscript denotes that dy, , , is the change 


in vapour composition which would occur if 


plate (N 
the column. 


1) were in the enriching section of 


Substituting for dy, ., from (2) and for dyy,, 
from the corresponding equation for the stripping 


section into (14) 
(1 — 


ap (1 


as | day +1 
E) ba vl 
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Transforming (15) and substituting (4) and (6) 
into the resulting equation 


4 
Dexp[(N 1) vy] 
(1 — E)[C exp (No,) 
D exp (Nov,)}} 


ag {Cexp[(N 


xp { A exp[(N | 
Bexp[(N + 1) ug] 
(1 A exp (Nu,) 


B exp (Nu,)}} (16) 


From the definition of v, , (7) and (12), it may be 
that the left 
terms proportional to 


wen 


seen hand side of (16) consists of 


exp 


and to 


exp (N 2M 1) / 


Ss 


Sp l Ss) 


In expressions of this type the Laplace operator 


p appears to have the same effect as an ordinary 
variable which lies along the >» Oas 
0=<T ©. Thus it 


initial part of the response, i.e. p apparently large, 


range © —p 
may be seen that for the 
the C terms in (16) are small compared with those 
D and 
from (9) it can be shown that the B terms in (16) 


containing may be neglected. Similarly 


are also negligible and the equation simplifies to 
agDiexp[(N 
xpAtexp|[(N 


1) | (1 


E) exp(Nov,)} 
E) exp(Nu,)} (17) 


Plate (N + 1) 


Similarly by considering a 


(iw) 
material balance 


over plate (N + 1) in the unsteady state 


LyAexp(Nu,)—Lg Dexp( Nov.) + FAry/p =0 (18) 


Solving for A and D from (17) and (18) and 
hence determining B and C from (9) and (11) and 
substituting these values into equations (4) and 
(6) gives 


12 
- = 2 
‘ 
1 + Ss 
+ 
15 
(15) 
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[exp (v,) (1 « 


(ay Ly | exp | 


ay pjl—(l—E)[1 


FAr, pil (1 (ag Ly ap Lg) exy 


(pwap Sy ag Sp) 4 [ E (xs 


exp 
(E (a5 —1)Ss 


(pwap Ss / ag Sp) 


Following WILKINSON and ARMSTRONG and usin 


ing equations are 


tp 
Na, 


exp | (N 4 1) ve] 


(‘ay Ly ‘ap Lg) exp — | 


1) Ss/as] + v, 


column to changes in feed composition 


‘xp | (N 


U,) — (4 apl.s) exp (v, — u,)} 


1) 


u, + (ap’ 1)/xp'] exp nu,) (19) 


xp (nu,) 


(— u,) 


»(— U,) — (ag Le ag Lg) exp (v, — u,)} 


exp[(M + —(M l (20) 


g the pseudo-equilibrium line [2] the correspond- 


| 
‘xp (nu,) Cap (21) 


Leg 
and 
FAa, — (‘ag Lp, ‘ap Ls) exp (v, uy) | 
nolL, IV Ve 
exp (mig) (prolp Ls) *8 Ls] -exp|(M 1) v, (MW m) (22) 


where 


LE 
E + (ap V/L,) (1 BE) 
E E 
and 


E + (agV/ Lg) (1 E) 


are the slopes of the pseudo-equilibrium lines and 
H. The expres 


sions for u,, are identical for both cases, but 


the non-dimensional time 
using the pseudo-equilibrium lines 


(So — 1) 
(Ss + 1) 


1\? 2Lep 
MN Ls (1 Ss) 


(contrast equation 7) 


when the plate efficiency is unity the two pairs 


of expressions are identical. 


In order to transform X,, and X,..,, back into 


1—(1 — E)[1 — Ls)]e 


the time domain it is necessary to invert functions 


ol the form 


1)? ag Sp | 
2(1 + Ss) ap Sel 


In order to simplify the inversion of X, and 


Xy.,, it was assumed that for the initial part 
S,)?/2(1 + Sg). This 


would appear to be valid for columns operating 


of the response p > (1 


at fairly high reflux ratios and has been justified 


elsewhere [3]. Hence u, , and v, , simplify to 


(Ss 1) 2 
se OE 


NV 


(Ss — 1) Ss p 
Ss +1 as Sp (1 + Ss) 


With this approximation the inversion of (19) 
and (20) becomes fairly straightforward and is 


effected by expanding 


xp (— u,) — (ag apy Lg) exp (v, — u,) 


| 
12 
| 
and 
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by the binomial theorem and inverting the ex- 
pression term by term. (All the relevant Laplace 


transforms may be found in CarsLaw and JAEGER 


[4}). 


Ext! NSION OF THE THEORY TO BIGGER VALUES 
or TIMI 


For the transient response of a single enriching 
section to a step change in feed composition, tt 
has been shown 15) that the solution is 


J, exp (23) 


Response of the top product composition to 


It is thought that for the complete column the a step change in feed composition 


form. For the Theoretical Witkinson) O run F.T.P.1, 


response would have a similar sete 
run F.T.P. 2. 
approach to equilibrium (23) may be written as 


dr, (T) 
exp p,7 


F(T) is the fractional approach to equilibrium. 
Hence from the calculated steady state change 


in plate composition the response may be 


extended to larger values of time. 


Fic. 4 Response of the composition of plate 8 to 


a step change in feed composition, 


Theoretical run F141 
O run F142. 


Fic. 2. Theoretical response of the top product and Fic. 5. Response of the bottom product composition 


plate Stoa step change in feed composition, toa step change in feed composition. 
—— Theoretical Wikinson O run F.B.P.1 


Using pseudo-equilibrium line for plates 0 and 8 
run F.B.P.2. 


° Using FE directly for plates 0 and 8 
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COMPARISON OF THEORY AND EXPERIMEN' 

The experimental results of Witkinson [6] 
were obtained at conditions which approximate 
closely to those required for comparison with the 
present theory. 

For the top product and plate 8, (19) was 
evaluated up to T 100. The theoretical results 
F,(T)| against T on 


semilog paper (Fig. 2). Good straight lines 


were then plotted as {1 


could be drawn through the points and therefore 
it was possible to extend the curves to larger 
values of time. The theoretical curves are plotted 
with the results of WILKINSON on Figs. 3 and 4. 
Similarly (21) derived using the pseudo-equilib- 
rium lines was evaluated for the top product and 
plate 8, and is also plotted on Figs. 3 and 4. 

As an example of the use of the equations for 
a plate in the lower half of the column, (20) was 
evaluated for the initial part of the response and 
may be compared with the experimental results 
in Fig. 5. 

Discussion 

The agreement between theory and experiment 
is within 20 per cent, which is thought to be 
reasonable when account is taken of the many 
approximations made in the theoretical derivation. 
The agreement between the curves obtained 
using the pseudo-equilibrium lines and plate 
efficiency directly is also within 20 per cent and 
thus perhaps the increased labour required to 
evaluate (19) compared with (21) is hardly justified 
under these operating conditions. At lower values 
of the plate efliciency the divergence is likely to 
be greater and the more accurate solution might 
then be warranted. 

The values obtained for the first decay constant 
p, for the top and bottom products and for plate 8 
are — -O071, 0083 and 


units of non-dimensional time. These values are 


0083 reciprocal 


quite close and give added justification for the 
method used to extend the solution as the column 
approaches equilibrium. 


CONCLUSIONS 
Theoretical expressions have been derived for 
the response of a distillation column to changes 
in feed composition. These expressions have been 


derived both by using the pseudo-equilibrium 


~ 
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lines and by using the plate efliciency directly. 
The theory agreed with the experimental results 
of WILKINSON to within 20 per cent which is also 
the maximum divergence between the two 


theoretical expressions, 


NOTATION 


x slope of equilibrium line 
x slope of pseudo-equilibrium line 


(1 — E)+E 
+ E(Lg 


 ELg + ag V(1 — EB) 
8 = constant in the equation of the vapour-—liquid 
equilibrium line y B 
m plate number stripping section 
n plate number enriching section 
Pp Laplace operator 
r a constant 
s a constant 
t time following a disturbance min 
t;,» = defined by (5) 
0,9 = defined by (7) 
u ratio of the reboiler hold-up to plate hold-up 
1 mole fraction of the more volatile component in 
the liquid 
y = mole fraction of the more volatile component in 


the vapour 
functions in equations (4) and (6) 


\ 
E = Murphree plate efficiency (liquid) 
I molar feed rate mole /min 
F(T) = fractional approach to equilibrium 
I] = plate liquid hold-up mole 
number of roots 
J = exponential multipliers 
L = liquid flow rate mole, min 
M number of plates in the stripping section 
N number of plates in the enriching section 
aV 
aV(l — E) + EL 
ELy + (1 BE) 
EH 


mole /min 


non-dimensional time 


V’ = vapour flow rate 
X = Laplace transform of 84 


Subscripts 
i = general term 


m = plate number stripping section 
n plate number enriching section 


E = enriching section 
F feed 
= stripping section 


~ 
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On the response of stagewise processes to stationary randomly- 
fluctuating inputs 


A. Acrivos 
Department of Chemical Engineering, University of California, Berkeley 
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Abstract—This theoretical paper deals with the response of stagewise processes to stationary 


randomly-fluctuating inputs with Gaussian probability distributions. It is shown how the problem 


may be formulated mathematically and solved in a manner believed to be more straightforward 


than the techniques presently used in control theory. Exact solutions are obtained for linear 


systems and two asymptotic series are developed for cases where the fluctuations are either 


slow or rapid. By means of these expressions, which involve only simple algebraic operations, it 
is possible to calculate without too much difficulty the probability that certain critical variables 


in the system will exceed a prespecified range of values. 


Résumé—Cet article de la réponse des systémes étagés A des signaux stationnaires préesentant 


des variations aléatoires suivant une loi de distribution gaussienne de la probabilité. Il montre 


comment le probléme peut étre exprimé sous forme mathématique et résolu de facon a obtenir 


une estimation de la réponse de maniére plus directe que par les méthodes actuellement utilisées 


en théorie du contrdle. Il indique les solutions exactes obtenues pour des systémes linéaires et 


donne deux développements cn séries asymtotique pour le cas de variations lentes ou rapides 


respectivement, Au moyen de ces expressions, qui ne comportent que de simples opérations 
il est possible de calculer sans grande difliculté la probabilité pour que certaines 


alyveébriques 
variables critiques du systéme dépassent des valeurs fixées 4 lavance. 


Zusammenfassung—Diese theoretische Arbeit behandelt die Antwortfunktion von stufenweisen 


Prozessen auf stationiire, zufillig schwankende FEingangsfunktionen mit Gaus schen 


Wahrscheinlichkeitsverteilungen. Es wird gezeigt, wie man das Problem mathematisch formulie- 


ren und wie man es zielstrebiger losen kann, als mit den gevenwiirtig benutzten Mitteln der 


Theorie der Regeltechnik. Fiir lineare Systeme werden exakte Ldsungen erhalten und fiir Fille, 


bei denen die Schwankungen entweder langsam oder schnell sind, werden zwei asymptotische 
Reihen entwickelt. Mit Hilfe dieser Ausdriicke, die nur einfache algebraische Operationen 
erfordern, ist es méglich, ohne zu grosse Schwierigkeiten die Wahrscheinlichkeit zu berechnen, 


mit der bestimmte kritische Variable im System einen vorgegebenen Bereich iiberschreiten. 


INTRODUCTION for example, it is known that in many polymeriza- 


THE response of stagewise processes to changes in _ tion reactions the reaction temperature may have 


the input variables has received an increasing a profound influence not only on the quality but 


amount of attention in the last few years, for it sometimes even on the nature of the product, 
is being realized that systems designed solely and it is for this reason that in such processes the 
on the basis of their steady state performance the temperature cannot be permitted to vary 
may possess undesirable transient characteristics more than a few degrees from its desired value. 

and may therefore be diflicult to control. In Such problems are of course quite common in 


particular, since it is usually highly advantageous many of the more important stagewise processes 
for the product of a given process to be of essen- in chemical engineering — stirred tank reactors, 


tially uniform quality, the operating system distillation columns, extraction units, etc. to 


including the controller must be so designed mention but a few-—and it is clear therefore 
that the characteristic properties of the products that before one can attach a suitable control 
should not be allowed to vary outside prespecified, system to an operating unit one must be in a 
position to predict how the process will respond 


and sometimes quite narrow, band limits. Thus, 


12 
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to input disturbances which are known or expected 
to occur. Such upsets are of two rather general 
types: they either take place in some orderly 


fashion and may therefore be represented by 


well-behaved mathematical functions of time, 
or they occur in a purely random manner which 
is describable only in a statistical sense. So far, 
the majority of papers on this subject which have 
appeared in chemical engineering literature 
(1, 2, 3} have studied transient phenome na of the 
first type with sinusoidal or step functions as 
input disturbances. In many practical situations, 


however, the input fluctuations are of a com- 
different 


unpredictable and random manner. 


pletely nature, for they occur in an 


The theoretical analysis to be pres nted below 
will deal precisely with this problem of how a 
stavewise S\ stem will be expected to rr spond to 


such randomly fluctuating input conditions. 


Mathe matically, the problem IS by no means new. 


Its importance was recognized some time ago 


by electrical engineers, and a commonly em 
ploved method of solution can be found in many 
of the standard texts on control theory [4, 5, 6}. 
On the hand, the of these 


established tools of analysis to stagewise systems 


other applicability 


of particular interest to chemical engineers was 
not demonstrated until quite rece ntly by Karz [7 
and by Arts and AmuNpbson [8]. It turns out, how- 
ever, that due to the considerabk complexity ot 
the typical chemical engineering stagewise 


the 


pro 


and therefore of mathematical model 


which describes its behaviour, the usual procedure 
for attacking this 


becomes quite cumbe rsome. 


rather involved problem 
The present study was therefore motivated by 
the 


deve lopment 


desire to devise a formal mathematical 


which would be somewhat more 
straightforward and easy to manipulate than the 
techniques presently used in control theory and, 
as a consequence, suitable for handling stagewise 


systems with a large number of components. 


T 


SYSTEMS 


Mop I 


ForMAI 


MATHEMATICA! LINEAR 


Its 


It has already been established [1, 8, 9] that, 


AND SOLUTION 


for small perturbations about the steady state, 


the transient behaviour of any stagewise system 
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is governed mathematically by the matrix equa- 
tion 

dx 

dt 


Az 2 (1) 


vector a 


the column }a;; has as its 


components the characteristic variables of the 


where 


measured from their 


state 


stagewise process corres- 


ponding steady values. In a system of 


stirred-tank reactors, for example, these variables 


would represent in general the temperatures, 
concentrations and flow-rate variations in the 
different interconnected units. On the other 
hand, the elements z,; of the column with z = { 2; | 


are directly related to the different input fluctua- 
tions which cause the system to depart from 
the desired steady state, while the components 
a,, of the constant matrix 4 


without difficulty 


may ns i 
the 
steady state solution of the appropriate algebraic 


rule be calculated from 
equations which characterize the process [9]. 


If now z were a known function of time, then 
the response of the system to such a disturbance 
could easily be obtained from the general solution 
of equation (1). If, however, as we shall suppose 
from now on, the elements of z, are all stationary 
random functions which may be described in a 
statistical sense only, then the problem at hand 
is indeed quite different in nature from the usual 


This 
theoretical 


transient phenomenon of control theory. 
the fact that 


analysis is supposed to yield in such a case is 


follows from what a 
not the instantaneous value of each component 
r,; in the system (which is clearly impossible since 
z cannot be completely specified in the usual 
sense) but the probability that one or more of 
the 


critical 


variables a 


; will lie outside predetermined 


limits. Thus, for example, if it were 
known that the reaction temperature in a given 
process should not be permitted to exceed a 
certain upper limit or fall below a specified value 
lest undesirable side products be formed, then it 
would indeed be of considerable interest to the 
control engineer to be able to predict the pro- 
bability of such an event occurring. 

To simplify the development we shall suppose 
next that all the elements z, of the column vector 


z are normally distributed with zero mean, as is 


: 

: 

: 
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usually the case to a first approximation at 
It then follows [5] that the elements of 


r are also Guassian. Therefore. the probability 


least. 


distribution of all the elements x,, and in particular 
the probability that they would lie outside set 
limits, could readily be calculated by the standard 
techniques of statistics if the variance of each 
r, were determined. Thus the main objective of 
our analysis will be to show how these variances 
may be derived from equation (1). 
Before proceeding with our solution however 
find it advantageous to introduce here the 
the ele- 


we 


square correlation matrix Z £7 


ments of which are defined by 


(h) = 2, (t) z(t + h) (2) 


he either with 


the taken 


respect to time, or, in view of the ergodic hypo- 


where average may 
thesis, in the probability space. We see then that 


2, (4) is the familiar cross-correlation function 


[4, 5] between the two random variables z, and 
z; which, in the absence of any periodicity in the 
h (4). This 


correlation matrix is a function of h only and 


fiuctuations, vanishes as > 00 
not of time, if, as we shall suppose, we are dealing 
only with stationary random disturbances, the 
mean properties of which in a statistical sense 
are time independent. We shall further suppose 
that Z (hk) is already known since it is related to 
the uncontrolled input upsets which cause the 
system to fluctuate. 

Let us next postmultiply each term in equation 
(1), at t, the 


vector Zp h) , Upon averaging we can 


which is evaluated time by row 
show that the modified equation (1), may be 


transformed into 


where : 


h) x, (t 


te 


At) 
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The solution to equation (3) subject to the 
boundary conditions 
p 0 for h + L (5) 
Is 
® — exp Ah | exp(— Ar) Z(r) dr (6) 
h 
This may be verified by direct substitution. 


Furthermore, since the process 1s assumed to be 
stable about its average steady state, the charac- 
teristic roots of 4A must have positive real parts 
{9} and thus, as required by equation (5), it is 
indeed true that 


lim exp Ah 0 


x 


On the other hand, if the transpose of equation 
(1) 

day 

dt 


which is evaluated at time (¢ h) is premultiplied 


through by 2 (t), and then averaged, it may be 
reduced into 
dX 
XA,=@ (8) 
dh 
where 
(t) a; (t h)} (9) 
will be recognized as the correlation matrix of 


the process output variables. Therefore, if again 
X Oas h—> ox. 
h 


| P(r) exp (Ay t)dr exp | 


x 


and, in view of equation (6), 


(t h) x, (t) 


At 
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+ AD = Z{(h) 
dh (3) 
(t) z(t + A), and 
(4) 


A. 


exp (Ar) 


. 


Finally, if. as is usually the case, one ts 


need oniv to compute 


r., then one 


equation (11) bv setting / 0. 


MATHEMATICAL SIMPLIFICATIO? 


ForMAI SOLUTION 


It has been established in the preceding section 


that the formal 
may be given by the rather impressive looking 
It can be ourse that, 
as expected, equation (12) is equivalent to the 
result that is obtained by the 


techniques, vet it is felt that not only ts our 
t that the 


solution of our linear problem 


equation (12). shown of 


more standard 


derivation more straightforward b 


form of the solution is more suitable for numeri al 
computations especially for compl ted systems 
Ihe reason for this 


products of 


with large order matrices. 
is that equation (12 which involves 


matrices only, is in general quite amenable to 


under 


manipulations especially 


mathematical 
those 
simplifications may be carried out 


conditions where, as we shall show, useful 


In general, the exact numerical evaluation of 


task. We 


portion ofl 


equation (12) Is a rather formidable 


shall therefore devote a considerable 


the remaining part of this paper to examining 


wavs in which we could simplify the rather 


involved mathematical expressions and thereby 
lighten any numerical work. To accomplish this 
in a straightforward manner we must first examune 
the integral 

x 

$ 


exp ( (13) 


Ay) Z (y) dy 


. 


note first 


which appears in equation (12). We 
of all that according to Sylvester's theorem [10)}, 
the matrix exponential function may be expanded 


in the form 
exp (- 
k=l 


exp (— Ay) Ay) U, (14) 


fy) Z (y) dy| exp (Ay 7) dr exp (— 


Ay) Z (y) dy 


(11) 


interested onlv in the variance of the individual fluctuations 


diagonal elements of the matrix Y,, which ts obtained from 


exp | 1,7) dr 


where n is the order of the matrix, A, its charac- 
teristic roots, assumed distinct and, in view of the 
postulated stability of our system, with positive 
real and U, are the idempotent matrices 


of A {10}, Therefore, 


pare 


. 


| exp | 


fu) Z (uy) dy 


eXp | \ Z (y) dy 


This equation is clearly of very limited usefulness 


for computational since it is” well 
known how difficult it is in general to evaluate 


the the 


matrices of A. It become 


puryp wes, 


numerically characteristic roots «wd 


idempotent will 
apparent, however, that equation (15) 1s indeed 


quite valuable, from a theoretical standpoint, 
since it will allow us to derive and to justify the 
approximations which follow. 

Let us then consider the typical element 2%, 
of the matrix Z. Since 2, (4) is the cross- 


correlation function of the two random variables 


square 


and 2, it has the general shape shown mn 
Fig. 1. In particular, we define the two positive 
characteristic parameters hi; and 0; by the 


following two relations : 


2, (A) « for all |h (16) 


and 
(17) 


where « is an arbitrary small number (c¢, for 
example, may be set equal to 0-1). Furthermore, 
if A, and A, are the characteristic roots of A with, 
respectively, the minimum and maximum real 


parts so that 
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Fic. 1. A plot of (a) typical cross-correlation function 


and (b) typical auto-correlation function of two 


random variables 


R(A,) = min R(A) 


(1S) 


R(A,) = max R(A,) 
j 
where the symbol R denotes the real part of a 
number, then we may segregate the random 
disturbances into two broad categories according 


to whether: 


l 
R(A,) > a, (19 
or 
R(A,) < h, (19a) 


We shall show now how equations (12) or (15) 
may be simplified if either of the above two 


asymptotic conditions is met. 


Case I 

If the input fluctuations are such that equation 
(19) is obeyed approximately by all the elements 
of the square matrix Z, then the integral given by 
equation (13) may be calculated conveniently 


as follows: We first expand Z (h) into a power 


series, so that 


of stagewise processes to stationary randomly-fluctuating inputs 


283 


Z (h) = Zy +hZ, + —Zy+..... 


where the representative kth term of the expan- 


sion is given by 


d* 


| ij 
| dh* 


It is then easy to verify that 


X, = AZ, Ay 
[4-*Z, Ap! Z, Ap] 
[ A-*Z, Ay! — A*Z, Ag 

A-'Z, Ap*)+... (20) 


0 


an expression which is indeed most convenient 
for numerical computations on a high speed 
computer, since it involves only products of 
matrices. Interestingly enough, though, equation 
(20) may also be derived by an iterative scheme 
based on our two fundamental equations (3) and 
(8). Thus, if to a first approximation we neglect 
both derivatives in these two equations, we find 


that 
@~ A X= A (21) 


and therefore 


X, = 4-1Z, Ag" 


which is identical to the first term of equation 
(20). If now the two derivatives are approxi- 
mated by making use of equation (21), and if 
the re sulting expressions are solved once more 
one can verify that indeed the first two terms of 
equation (20) are obtained. This iterative pro- 
cedure may then be repeated successively to 
vield equation (20) 1m its entirety. 

We see, therefore, that since the derivative 
terms in equations (3) and (8) represent in a sense 
the capacity of the stagewise system and are, 
as a consequence, related to its ability to dampen 
out the input fluctuations, the solution given by 
equation (20) is nothing more than an expansion 
about the so-called “* quasi steady state.” That is 
we first assume that the system can respond 
instantaneously to the applied disturbance and 
then correct our mathematical solution by an 
iterative procedure of successive approximations. 
Thus we would expect that this series expansion 


given by equation (20) would converge rapidly, 


j 
3 
| 
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from a computational standpoint, if the dominant 
of the input fluctuations were small compared 
to the equivalent residence time of the stagewise 
system. It can be shown that this is mathe- 
matically equivalent to the inequality given by 
equation (19), and we shall refer to such distur- 


bances as “ slow.” 


Case Ul 


On the other hand, if the dominant frequencies 
of the input fluctuations are large compared to 
the equivalent residence time of the system, and 


if in particular the input disturbances approach 


a mathematical idealization known as “ white 
noise 4, 8]. then the “ quasi steady state 
approximation is clearly unsatisfactory. The 


series expansion for Re given by equation (20) 
will, therefore, either converge very slowly or will 
even diverge, so that a different asymptotic 
solution is needed, 

We can now simplify equation (12) by noting 
that when all the fluctuations are of the above 
type, the various cross-correlation functions 

(h) approach zero even for small values of h. 
Consequently, since h,, - 0 (see Fig. 1), equation 
(19a) must be satisfied. This then suggests the 
following expansion of our solution. If we first 


rearrange equation (12) into 


. 


| exp (dr) Lexp (A, r)dr 
o 


[exp | | exp | Ay) Z (y) dy| exp (Ay dr 


. 


we can show that by expanding the exponential 
functions in the second term, 
o 


where the matrix L is given by 


L exp Ar) Z(r)dr 
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It is easy to verify now that indeed the conver- 
gence of equation (21) is very rapid if hy +0 
for all the elements 2, (h) of the matrix Z. 
Furthermore, L_ itself may be conveniently 
calculated by the series 


L Zir)dr A | 


x 


ePZir)dr+.... (23) 


. 
x 


- 


which again will converge quickly if all hy - 0, 

Unfortunately, however, the evaluation of the 
first term in equation (21) is somewhat involved 
and must be carried out by an indirect method, 
It is clear first of all that if the input disturbances 
belong truly to the type known as “ white 


noise so that 


(h) = 1, 


i 


where 4(h) is the familiar delta function of 
mathematical physics, the equation (21) would 
reduce to 


Xy exp (dr) Lexp Ayr) dr 


. 
x 


We see, therefore, that the first term in equation 
(21) represents nothing more than the “ white 
noise ’’ approximation to the input fluctuations, 


In that case, equation (6) may be simplified into 
® — exp(Ah)L for h <0 


and as can be verified by direct substitution into 
equation (8), 


X = exp (Ah) X*, for h = 0 


where the constant matrix X*, satisfies the 
algebraic expression 


AX*, + X% Ae =L | (24) 


The following is a rather straightforward scheme 
for obtaining X*,: If the transpose of equation 


(24) is added to equation (24) one can show that 


A[X*, + (X*%)7] 


[X*, A, L 7 L T (26) 


; 
|| 
= 
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We next let 


A[X*, + (X*,)7] =(L+ Ly) + F (27) 


where F is an unknown skew-symmetric matrix 


of the form 


hie Sis - | 
his fog 


Therefore, 
X*, (X*,)4 Ly) 1-1 F (29) 


The elements of F are finally computed by 
taking into account the fact that X*, + (X*,), 


is symmetric. Thus, by equating the approximate 


terms in equation (29) we can obtain aw L)/2] 


algebraic equations for the unknowns f/f... It 
should also be noted carefully that since, for 
reasons explained earlier, one will in general b 
interested only in the diagonal elements of x¢.. 
which are identical to the variances of the output 
variables, one need not proceed any further once 
equation (26) has been solved. 

It is clear then that equation (21) may be 


written in the form 


+4 | (r) dz 


where Xx, is given by equation (24). rhe com 
putation of this expression is quite straight 
forward once X*, has been evaluated, for the 
remaining terms in the series involved only 
products of matrices. Furthermore, since Z 
may usually be approximated by a linear com- 
bination of simple matrix functions, the integra- 
tions in equations (22) and (30) may as a rule be 


carried out analytically without much difliculty. 


The Lene ral case 


In any actual situation, the input fluctuations 


willnaturally be of both kinds ; “slow” and “rapid 


In view of the linearity of our basic equations, 


however, the solution of this more realistic 
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problem may easily be derived. Thus if we let 

Z=2,+ 2, 
where all the elements of Z, and Z. satisfy. 
approximately, equations (19) and (19a) re- 


spectiy ely, then 


where X, is given by equation (20) with Z Z, 
and AX, is given by equation (30) with Z } 


CONCLUSIONS 


In this theoretical analysis we have studied 
the response of linear stage w ist systems to 


stationary randomly fluctuating inputs with 
Gaussian probability distributions. formal 
mathematical development of the problem was 
presented which is somewhat more straight- 
forward and easy to manipulate than the tech- 
niques presently used in control theory, since it 
leads to a solution which may be drastically 


simplified under asymptotic conditions. Thus an 


expansion about the ~ quasi steady state sis 
results in equation (20). whereas for ; rapid is 
iluctuations equation (30) is preferable. In 


general, the essential features of such a transient 
problem may be clearly brought out by segregat- 
ing the disturbances into two broad classes. 
‘slow ” and “rapid.” Then, using the appro- 
priate asymptotic expansion, which were derived 
in this paper for these two types of fluctuations, 
the computation of the variances of the critical 
output variables may be caried out without 
much difficulty, for it involves only simple 
algebraic operations suitable for a high speed 


digital computer. 


APPENDIX 
1 numerical example 


A simple numerical example will serve to illustrate the 
essence of our theoretical analvsis We shall consider a 
system which may be represented mathematically by 
matrix equations which involve only second order matrices 
for, as already shown by Arts and AmuNpDson [8], the 
transient response of a single stirred-tank reactor may be 


adequately described by precisely such a mathematical 


model, 
Let 
2 1/3 1 O 
1 | and Z exp ( a ih) | 
l 1/2 0 0 


| | 
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\. AcrIVOs 


where « is a positive parameter. By performing the 


appropriate calculations we can show that 


O-3035 2-19065 


0- 1038 O-1761 O-1761 


without difficulty 
finds that 


integrated 
Cone 


Equation (12) can now be 


by the use of Sylvester's theorem 


| 


O-O10T7 


OA, 0-05484 


4187 


O473: 


(A, O47355 O 2791 


This exact solution has two asymptotic forms dep« nding 
on whether «— 0 or 
Case I. The 


equation (20), is applicable when « —>- 0 


state approximation 


Thus by 


* quasi steady 


neglect 


ing all but the first term in equation (20), we can show that 


0 16 
(b) 


must be considered when 


Case 
In particular if only the first tern 


Equation (30) 


o of equations 
(23) and (30) are taken into account, then 


O50 


0-300 ooo) 


a 
Fic. 2 


function of « 


as 


The variance of the output variable 2, 


for the numerical example 


2,, (0) and a,, (0), respectively the variances of the two 
random output variables ry and ty, are shown in Figs. 2 
and 3 as functions of the parameter a, together with their 
and (c). It 
is clearly seen that both ry, (0) and oq (0) may be 
quately described by the leading terms of their asymptotic 
It is felt, therefore, 


asymptotic forms as given by equations (b) 


ade- 


expansions about « O and «a— @, 


a 


3. The variance of the output variabl tr usa 


function of « for the numerical example 


that this simple numerical example has served to illustrate 
our main point of view, namely that for those high order 


processes for which an exact solution cannot be con 


veniently obtained, the two asymptotic expansions given 
by equations (20) and (30) will indeed provide in general 
an acceptable representation for the main characteristic 


features of the system. 


NOTATION 


the independent variable in the cross-correlation 


functions (see equation 2) 


a;; |}, & constant square matrix, characteristic of 


the system (see equation 1) 
r;}, a column matrix the components of which 


are the characteristic variables of the system 


z;}, a column matrix with the fluctuating input 


variables as its elements 


(t+ z; (0), the cross-correlation function of 


the two random variables z; and z, 
=, + (0), the correlation matrix of the 
inputs 


(tat A)}, the correlation matrix of the 


variables of the system 


dj; (h)} where ry (t) (t h) 


/;;}, a square matrix defined by equation (22 
the characteristic roots of A {10} 
the idempotent matrices of A [10] 
the order of the matrix A 
characteristic parameters of the cross-correlation 
function zi (h) (see Fig. 1) 
Subscripts 
0 = denotes the matrix evaluated at h = 0 


T = the transpose of a matrix or a column vector 


j | 
“oO 
1 2a + Ay + Ay — 0 
Ag (A, a) (A, x) 5582 
(a) 
vol 
~) L2 
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Influence of mass transfer on coalescence of drops 


Received 29 January 1960) 


IN RECENT years increasing interest has been shown der 


interfacial phemonena during transfer between 


phases. Lewis {1} found that at liquid-liquid interfaces 
severe disturbances rhe ex- 
that 


lients im the 


solute transfer may cause 
planation offered by Stawart and Nassensrein [2 
these phenomena occur as a result of gr 


interfacial tension caused by non-uniform iss transfer 


along the interface seems to be satisfactory 


The above studies have been made by observations of 
flat interfaces or single drops suspended in a continuous 
liquid phase. Less attention has been paid to the in- 
dustrially more important case of the behaviour of swarms 
JOUNSON 


that in a 


of droplets under the influence of mass transfer 


and Buiss [3] as early as 1946 found spray 
when the 
to the 
When solute was extracted from the drops much 
that 


drops. Pratl 


extraction column transfer was highest 


THLASS 


solute was transferred from the continuous 


phase 
droplets 


larger drop sizes were observed, a phenomenon was 


ribed to ime reased coule sciny of and 


co-workers [4] found that in different types of extraction 


equipment, droplets considerably larger if the 


were 


solute was transferred from the droplets than if transfer 


was in the opposite direction As a result, maximum 


allowable flow rates in these cases were very hig For 


example, with toluene as the dispersed phase the flooding 
rate of a rotating disk contactor column is 100-200 per cent 
higher when acetone is transferred from the droplets to 
the continuous water phase than when the direction of 
transfer is reversed or when no soluté 


The 


ascribed to droplet oscillation as found in the single drop 


is present at all [4 


increased coalescence rate has sometimes been 


studies mentioned earlier. It has however been shown in 


a previous communication that the actual cause of droplet 
oscillation, viz. the development of interfacial tension 
gradients is also responsible for the differences in stability 
of foams in the related field of gas-liquid contacting [5]. 


The explanation given may be repeated here shortly 


when during a mass transfer operation two gas bubbles 


in a liquid approach each other, the small amount of liquid 
between the bubbles tends to come into equilibrium with 


of the 
Phe compound 


the gas phase very quickly, while along the rest 
bubble area mass transfer is still going on. 
diffusing from the bubble phase will therefore be present 
at the highest concentration in the liquid film between the 
bubbles. When this higher concentration is accompanied 
by a lower surface tension, liquid is drawn away from the 
interfacial film effect) 


In the opposite case ¢ oalescence is retarded and 


(Maragnoni and coalescence is 
promoted 


the gas bubbles are very stable (foaming) 


occurrence of coalescence in liquid-liquid systems. 


It appears this explanation also holds for the 
This 
may be illustrated by considering the system benzene—acet ic 
this 


We assume for example that 


acid-water. Approximate interfacial tensions for 


system are shown in Fig. 1. 
the benzene phase 


is dispersed, and that it originally 


contains 1 per cent of acetic acid. According to Fig. 1, 


warren 


equilibrium compositions ( 


Approximate interfacial tensions at 20 °C for 
vol.) in the system : 


benzene—acetic acid-water (dyn/em.) 


the interfacial tension with respect to the small amount 
of equilibrated water phase in the contact area between 
of about 
Along the rest of the drop, mass transfer 


adjacent drops will rapidly approach a value 
20 dvn 
continues, the relatively small amount of liquid in the 
drop now tending to come into equilibrium with the large 
amount of surrounding water phase, which is solute-free. 
Interfacial tension will therefore continuously increase 
and finally approach the value of 35 dyn/cm. for the 
solute-free system. Under the influence of the different 
surface forces, water will be withdrawn from the contact 
area, thus promoting coalescence. 

It is easily shown with the aid of Fig. 1 that when water 
is the dispersed phase and acetic acid is transferred from 
the drops to a continuous benzene phase, the interfacial 


tension gradients developing around the drops have the 


28S 


ae 

bs 

7 

ACETIC & 

\ 12 
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‘ 

4 

é 


Bonzem acid water system acid 


transfer from drops to continuous water phase immediate coalescence 


hia. jonzene Cl, acetic acid-water svstem. Acetic acid transfer 


from continuous water phase to drops, strong retardation of coalescence 


an 
| 


vol 
12 
4 

: 

: 


same sign and coalescence is also promoted. Upon the ad- 
dition to two immiscible liquids of increasing amounts 
of a third liquid, miscible with both, miscibility will in 
general be increased. Therefore the tendency of interfacial 
tension change with solute concentration is generally the 
same as that shown in Fig. 1. The statement [4) that 
coalescence occurs whenever there is mass transfer from 
the drops into the continuous phase is therefore in accor- 
dance with the explanation offered here 

When, on the other hand, solute is transferred to the 
droplet phase, the above reasoning predicts that coales- 
cence will be retarded and the emulsion will become 
relatively stable. This is again illustrated by reference to 
Fig. 1. When the aqueous solution around the benzene 
drops contains e.g. 5 per cent acetic acid, the interfacial 
tension in the whole system will eventually approach about 
24dyn_ em (equilibration of a relatively small droplet with 
an excess of continuous phase), but the first thing to happen 
will be an advanced extraction of the water film in the 
contact area between drops (equilibration of a relatively 
small volume of continuous phase with an excess of droplet 
phase); this will increase the interfacial tension locally to 
about 35 dyn/cm. Therefore, continuous phase liquid will 
be drawn towards the contact area, stabilizing the drops 

Since little evidence of retardation of drop coalescence 
can be found in literature, experiments were carried out 
to contirm the above point of view. In these experiments 
drops were caused to form at the ends of two capillary 
tubes. These tube ends were facing each other in a beaker 
filled with the liquid which served as the continuous phase 
During their growth the drops approached each other 


The time passing between the moment of first contact and 


the occurrence of coalescence was observed 


Lewis J. B. Chem 


2) Srowarr and Nassensrerm H 


Jounson H and Biiss H 


Zuiperwec F. J. and Harwens A 


PRAT! 


Influence of mass transfer on coalescence of drops 


Rer 


Enong. Sci 


Trans 


R. C. e al Trans. Inst. Chem. Engrs. (Lond). 1957 35 301 


The following ternary systems were used : 
benzene—acetic acid-—water ; 

carbon tetrachloride—acetic acid—water : 

carbon tetrachloride /benzene—acetic acid—water : 
methyl isobutyl ketone—acetic acid—water : 


chloroform—acetone water. 


In all systems the liquid mentioned secondly was used 
as the solute. 

Without a solute the time needed for coalescence 
varied widely within each of the systems investigated. 
The mean values for the various systems were found to 
show marked differences as well. With a solute. however. 
the behaviour was quite different. When the solute was 
dissolved in the drops, coalescence invariably occurred 
almost at once after the first visible contact. irrespective 
f whether the organic solvent or the water was the dis- 
persed phase. With the solute dissolved in the continuous 
phase, however, coalescence was strongly retarded and 
in some cases did not occur at all, not even if big drops 
had a large area of contact (see the photographs in Fig. 2). 
The same behaviour was observed when both phases 
contained a solute, with an excess (above the equilibrium 
value) in one of them. The experiments thus completely 
contirm the analogy between the effect of mass transfer 
on dispersed phase interaction in gas-liquid systems and 
that in liquid-liquid systems. The authors hope to report 


on more quantitative aspects of the latter case in a future 


paper 
Koninklijke, Shell-Laboratorium, H. Grooruvuts 
fmsterdam J. Zomerwes 
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Edited 
256 pp. 


Techniques of Polymer Characterization. 
by P. W. 
Os. 


Auten. Butterworths, London 1959 


Tuts book presents a survey of eight major laboratory 
techniques used in polymer science in as many chapters 
written by academic and industrial experts. P. W. ALLEN 
writes briefly about choice of solvents and preparation of 
solutions ; R. W. Haut 


on osmometry ; D. 1 


on fractionation ; H. T. Hookway 


RusHMAN on other colligative 
molecular weight measurements ; F. W. Peaker on light 
OnYON on viscometry ; G. F. Prick on 


and R. J 


of graft and block polymne rs 


scattering ; P. F 
haracterization 
wndard of the 


end-group analysis ; CERESA on « 


The general st 


book ts acceptabk and occasionally comes near to per- 


fection, especially in its presentation of the theory and 


practice of light scattering by Dr. Prakes who has 


incorporated much of his own unrivalled experience in a 
Several hundred recent 


well-balanced historical account 


literature references make the book as a whole a valuable 
source of technical information 

Some rather minor sources of irritation must be men- 
tioned. In « book of this kind ther 
ences to materials, their trade-names and suppliers 


they 


must be many refer- 
Here 
are given with a disconcerting lack of uniform clarity 
Araldite 
mended without any such details, while 
Membrantilter Gesellschaft 


and completene ss. resin, for instance, is recom- 


i Bakelite varnish 
is merely referred to as V9700 
101 is amplified to Membranfilter 
Federal 


Institute 


(Gottingen) (sic) on p 
Gesellschaft 
p- 103. The address of the 


in connexion with 


(Gottingen, German Republic) on 


Shirley is not given, 


lest Leaflet 


* (mentioned on p. 230). The fullest possible 


though equally relevant 
No. Chem, 6 
documentation on anvthing that is worth recommending 
or quoting would surely be appreciated by readers of any 


book concerned with techniques. Incidentally, the purifica- 


Book 


Proceedings of the 1960 Heat Transfer 


and W. G. Vincenti. Stanford University Press 1960 


und Fluid Mechanics Institute 


tion of reagent-grade benzene (p. 125) can now be avoided 
by buying the grade suitable for mole« ular weight measure- 
ments (m.p. 5-4 °C; from British Drug Houses, Poole, 
Dorset). 

The description of transport of solvent through a semi- 
is to be 


permeable membrane as * liquid flow’ (p. 91) 


deprecated ; so is reference to the term c* as the cubic 


term (p. 109). The reader is warned (p. 82) that even preci- 
sion-bore glass tubing may vary in diameter by 0-1 mm 
along a single length of capillary. The reviewer has never 
this, though the 


limit 


found variations nearly as large as 


manufacturers. in their commendable caution, 


their guarantee to such tolerances. In practice one finds 


that such tubing is often as good as that des ribed by 


F. S. Darvron et al. (Trans. Faraday Soc. 1957 53 1269) 


which varied only by 00005 mm over a length of 


One effect of reading this book is to feel one’s appetite 
Orie must be felt 


whetted for disappoimtment 


that 


hore 
a more comprehensive and ambitious book has not 


emerged. Over 10 years ago, Weisspercer’s Volume I 


of “ Physical Methods of Organic Chemistry,” though not 


restricted to science, did not include 
admirable chapters (now out of date) on viscometry and 


but besides devoted 110 pages to ultra- 


polymer only 


osmotic pressure 
centrifuge techniques. Surely Dr. ALLEN’s pleas of lack 
of space, and of lack of general use of the ultracentrifuge 
will be widely regretted and challenged. Let us be grateful 
that eight experts can be found within 100 miles of Oxford 
to bring up to date the exciting story of polymer character- 
fields 

‘even in its present form, will save many 
But let us 


hope for a much fuller second edition in due course. 


ization in the included. “ Techniques of Polymer 
Characterization 


man-hours at the polymer scientist’s bench, 


N. Gorpon. 
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